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Q . Abstract: We systematically explore the spectrum of gravitational perturbations in 

^ . codimension-1 DGP braneworlds, and find a AD ghost on the self-accelerating branch 

of solutions. The ghost appears for any value of the brane tension, although depending 
on the sign of the tension it is either the helicity-0 component of the lightest local- 
Oh. ized massive tensor of mass < m? < 2H^ for positive tension, the scalar 'radion' 

^ , for negative tension, or their admixture for vanishing tension. Because the ghost is 

gravitationally coupled to the brane-localized matter, the self-accelerating solutions 

KA \ are not a reliable benchmark for cosmic acceleration driven by gravity modified in the 

IR. In contrast, the normal branch of solutions is ghost-free, and so these solutions 
are perturbatively safe at large distance scales. We further find that when the Z2 orb- 
ifold symmetry is broken, new tachyonic instabilities, which are much milder than the 
ghosts, appear on the self-accelerating branch. Finally, using exact gravitational shock 
waves we analyze what happens if we relax boundary conditions at infinity. We find 
that non-normalizable bulk modes, if interpreted as AD phenomena, may open the door 
to new ghost-like excitations. 
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1. Introduction 

Ever since Einstein introduced his famous "biggest blunder", the cosmological constant 
has been one of the most frustrating, yet intriguing aspects of General Relativity (GR). 
Ironically, just as Einstein needed a A to make a static universe, if we take his theory 
of GR as the description of gravity at the largest scales, we now seem to need a A to 
account for the cosmic acceleration observed at redshifts z < 1.7 [1, 2, 3]. Unfortunately, 
manufacturing a sufficiently small, positive cosmological constant from a consistent 
theory is not entirely straightforward, to say the least. The methods of effective field 
theory have so far failed to yield a satisfactory microscopic theory of the cosmological 
constant [4, 5]. Moreover, while the mystery of the cosmological constant is usually 
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posed as a problem for the field theory of matter, one may even wonder if in fact 
it might really be related to our formulation of gravity and inertia. Our hands-on 
experimental knowledge of gravity conforms with GR at distances between ~ 0.1 mm 
[6, 7] and, say, ~ 10 — 100 MPc. At these large scales we enter the domain of dark 
matter, a necessary component of the standard cosmological model needed to explain 
galactic rotation curves, which cannot be accounted for with GR and baryonic matter 
alone. At the moment, dark matter still needs to be completely explained by particle 
physics despite a plethora of reasonable candidates. A popular common theme in recent 
research is that perhaps it is not matter that is needed, but a modification of Newton's 
law and/or gravity at large scales. This idea is not new: ever since galactic rotation 
curves were found to be inconsistent with the luminous matter, such alternatives have 
been pursued [8]. 

While it is natural to hope that modifying gravity could be an interesting alter- 
native to dark matter, why might one hope that it could help with the cosmological 
constant? To illustrate this, we offer the following simple, heuristic argument. It is 
clear that in the Einstein-Hilbert action, the cosmological constant term appears as the 
Legendre transform of the field variable \/—g = a/| (iet(g^ 



JflU J 



From the canonical field theory rules, this means that this term trades the independent 
field variable A/|det(^^^,^ for another independent variable A. This is exactly the same 
as in quantum field theory, where one defines the generating function of the theory 
by shifting the Lagrangian by a 'coupling' J (j){x)J{x). This trades the independent 
variable for another independent variable J. After this transformation, the variable 
J is not calculable; it is an external parameter that must be fixed by hand at the end of 
the calculation, by a choice of boundary conditions. Once J is fixed to some value, is 
calculable in terms of it. The only difference between the usual field theory Legendre 
transform and the cosmological constant term arises because of gauge symmetries of 
GR, which render ^^/\det{g^JJ^\ non-propagating. It is a pure gauge variable that can 
always be set to a constant number by a change of coordinates. Therefore the Legendre 
transformation (1.1) loses information about only one number, which must be fixed 
externally: namely, by the value of A itself. As a result, in GR the cosmological 
constant is a boundary condition rather than a calculable quantity. One may then hope 
that by changing gravitational dynamics one could render ^y\ det{g^^)\ propagating, so 
that, in turn, A is also rendered dynamical. This could provide us with new avenues 
for relaxing the value of A. Such hopes have been already expressed before on a few 
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occasions [9, 10, 11]. However, analyzing modifications of gravity systematically, to 
check if they remain compliant with the tests of GR, hasn't been easy. 

On the other hand, in recent years the braneworld paradigm has emerged as a com- 
pelling alternative to standard Kaluza-Klein (KK) methods of hiding extra dimensions 
and a new framework for solving the hierarchy problem. In this approach our uni- 
verse is realized as a slice, or submanifold, of a higher dimensional spacetime. Unlike 
in KK compactifications, where the extra dimensions are small and compact, in the 
braneworld approach they can be relatively large [12, 13], or infinite [14, 15]. We do 
not directly see them since we are confined to our braneworld, rather, their presence is 
felt via corrections to Newton's law. Many of the more fascinating phenomenological 
features of these braneworld scenarios arise in models of warped compactification. In 
warped compactifications the scale factor of a four- dimensional brane universe actu- 
ally varies throughout the extra dimensions, providing us with a new way of making a 
higher dimensional world appear four-dimensional. In general, one can conceal extra 
dimensions from low energy probes by either 1) making the degrees of freedom which 
propagate through the extra dimensions very massive so as to cut the corrections to 
Newton's law off at long distances, or 2) suppressing the couplings of the higher di- 
mensional modes to ordinary matter so that the 4D gravitational couplings dominate, 
ensuring that the corrections to Newton's law are very small at long distances. The 
latter case is naturally realized in warped models, so that even infinite extra dimensions 
may be hidden to currently available probes. 

Braneworlds provide a natural relativistic framework for exploring means of mod- 
ifying gravity. It was quickly realized that by using free negative tension branes, one 
could alter Newton's constant at large scales [16]. More dramatically, Gregory, Rubakov 
and Sibiryakov (GRS) [17] noticed that by combining negative tension branes with in- 
finite extra dimensions, it was possible to "open-up" extra dimensions at very large 
scales, making gravity effectively higher-dimensional very far away. However, it was 
soon discovered by the authors that this model contained ghosts [18]. This was unfor- 
tunate since the metastable graviton had many desirable gravitational properties, but 
from a particle physics point of view the existence of a ghost is disastrous. Soon after, a 
radically new braneworld model was put forward, the DGP (Dvali-Gabadadze-Porrati) 
model [19], with graviton kinetic terms on the brane as well as in the bulk. The simpler 
versions of this theory are described by the action 

-^ 4 + n L' bulk Jbrane 



'{A + n)D bulk 

+ extrinsic curvature terms -|- / /^matter- (1-2) 

'brane 
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In general, there may be additional terms in the bulk. The key new ingredient here 
is the induced curvature on the brane. It could be generated, as it was claimed ini- 
tially, by quantum corrections from matter loops on the brane [20], or again in a 
purely classical picture of a finite width domain wall^ as corrections to the pure tension 
Dirac-Nambu-Goto brane action [22, 23]. Furthermore, it is also intriguing to note 
that induced curvature terms appear quite generically in junction conditions of higher 
codimension branes when considering natural generalizations of Einstein gravity [24] 
as well as in string theory compactifications [25]. Using holographic renormalization 
group arguments [26], DGP was shown to be equivalent in the infrared to GRS, how- 
ever, crucially, it appeared to be ghost-free, corroborating the perturbative analysis of 
[27]. This made it seem a real candidate for a new gravitational phenomenology at 
large distances. The induced curvature term yields a particularly interesting new phe- 
nomenon. In the case of a brane in 5D Minkowski bulk it allows for a self- accelerating 
cosmological solution [28], for which the vacuum brane is de Sitter space, with constant 
Hubble parameter H = 2M^/Ml, even though the brane tension vanishes. 

Clearly, the possibility of a fully consistent explanation of large scale acceleration 
is extremely exciting. It has generated a great deal of activity and investigation into 
the DGP set-up [29] (for a recent review, see [30]), with an astrophysical emphasis on 
black hole solutions [31], solar system tests [32], shock wave limits [33], and of course, 
whether DGP can truly explain dark energy [34]. Although many cosmologists have 
already embraced the DGP model, it has been found to suffer from various problems. 
There is the issue of strong coupling [35, 36, 37, 38], related to the feature that the 
graviton interactions go nonlinear at intermediate scales. More importantly, various 
investigations pointed out that there are ghosts on the self-accelerating branch [37, 38, 
39], however this debate still persists. 

Our aim here is to explore this issue in full detail. Since most of the explicit 
work on DGP has been done for the simplest case of a brane in fiat 5D bulk, with 
the dynamics given by (1.2), we will work in the same environment, and start with 
a review of this case. We will next consider the spectrum of small perturbations of 
the cosmological vacua of DGP, which describe a 4D de Sitter geometry. One of the 
confusing aspects of the literature on these braneworld perturbations (as opposed to 
braneworld cosmological perturbations!) is the alternate approaches of direct 'hands- 
on' calculations, which analyze the curved space wave operator for the gravitational 
perturbation directly [14, 40, 41] and the "effective action" approach, which was used 
to particular effect to confirm the ghost [18] of the GRS model via a radion mode 
analysis [42]. Naturally these approaches should be entirely equivalent, and we will 



^Harking back to the early manifestations of brancworlds [21]. 
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indeed see that. The technical comphcations in the identification of the spectrum of 
DGP gravity arise from the mixed boundary conditions for perturbations that may 
obscure the computation of the norm. 

The relevant modes in the spectrum of perturbations for addressing the concerns 
about stability are the tensors and the scalars. By going to a unitary gauge, we will 
see that the tensors are generically organized as a gapped continuum of transverse- 
traceless tensor modes, with 5 polarizations per mass level, and an isolated localized 
normalizable tensor, which lies below the gap. On the normal branch, this localized 
tensor is massless, implying that it has only two helicity-2 polarizations; on the self- 
accelerating branch it is massive, with < m'^ < 2H^ for positive tension, and has 
5 polarizations. When the brane tension is positive, the helicity-0 mode is the ghost, 
precisely because its mass sits in the region prohibited by unitarity, explored in [43, 
44, 45, 46]. Furthermore, the propagating scalar mode, or the 'radion', is tachyonic. 
This tachyonic instability of scalar perturbations is very generic, and by and large 
benign (see section 3.2). Moreover, the tachyonic scalar completely decouples on the 
normal branch in the limiting case where the bulk ends on the horizon^. On the 
self-accelerating branch, the scalar mode remains tachyonic but mostly harmless for 
positive tension branes, but as the tension vanishes it mixes with the helicity-0 tensor, 
and prevents the ghost from decoupling even in the vacuum by breaking the accidental 
symmetry of the massive tensor theory in de Sitter space in the limit m'^ = 2H^, studied 
in partially massless theories [43, 45, 46]. This mode becomes a pure and unadulterated 
ghost when the brane tension is negative, because it contains the brane Goldstone mode 
which does not decouple in a way similar to the GRS model [17, 18], consistent with the 
claim of [37]. Thus the self- accelerating solutions always have a ghost, and therefore do 
not represent a reliable benchmark for an accelerating universe in their present form. 
On the other hand the normal solutions are ghost-free, and thus may be useful as a 
model of gravitational modifications during cosmic acceleration. 

Our calculations further allow us to extend the analysis to perturbations which are 
not Z2-symmetric around DGP branes. This symmetry can be relaxed for braneworld 
models^, and actually this may be a more natural setting for the DGP setup which 
is more closely analogous to finite width defects or quantum corrected walls. In 
fact, in general braneworld models, when the requirement of Z2-symmetry is dropped 
one can get a whole range of interesting gravitational phenomenology, including self- 
acceleration, without appealing to induced gravity [48, 49]. We will show here that 



^Wc remind the reader that the situation here is similar to the single positive tension brane in the 
RS model where the radion also decouples. 

■^The Randall-Sundrum model [13, 14] was Z2-symmetric by construction, enabling to interpret it 
as a dual AdS/CFT with a UV cut-ofF and coupling to gravity [47]. 
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if Z2-antisyninietric modes are allowed, then in addition to the ghost, there is an ex- 
tra excitation which corresponds to the free motion of the DGP brane. This mode is 
tachyonic, and while it decouples on the normal branch in the single brane limit, it 
persists on the self-accelerating branch of DGP solutions. Nevertheless it still remains 
tame, since the scale of instability is controlled by the Hubble parameter, and so the 
instability may remain very slow. 

The presence of the ghost in the 4D description of the self-accelerating solutions of 
DGP indicates that the instability originates from the 'reduction' of the theory, and may 
not really represent a fundamental problem of the bulk set-up. A different prescription 
for boundary conditions might be able to circumvent the contributions from the brane 
localized ghost. However this requires rather special boundary conditions very far from 
the brane mass, that would not normally arise dynamically in a local theory. They 
allow the leakage of energy to, or from, infinity. Worse yet, an explicit exploration of 
potentials of relativistic sources shows that in this case other modes behave like ghosts, 
if interpreted in the 4D language. We can see this directly from the gravitational shock 
wave solutions which include the contributions from the modes that are not localized 
on the brane. 

The paper is organized as follows. In the next section we will review some of the 
salient features of the DGP model, describing its two branches of background solutions, 
the normal branch and the self- accelerating branch. In section 3, we will discuss the 
perturbation theory around the AD cosmological vacua of DGP, and identify its occult 
sector by an explicit calculation. In section 4, using gravitational shock waves, we will 
consider what happens when we include the contributions from non-normalizable bulk 
modes to the long range gravitational potential of brane masses. We will summarize in 
section 5. 

2. What are DGP braneworlds? 

We will work with the simplest and most explicit incarnation of DGP, where our uni- 
verse is a single 3-brane embedded in a 5 dimensional bulk spacetime. The bulk is 
locally Minkowski and the brane carries the curvature of the induced metric as well 
as the brane localized matter. The induced curvature terms will generically arise from 
the finite brane width corrections. The brane may be viewed as a (5-function source 
in the bulk Einstein equations, whose dynamics ensues from the total stress-energy 
conservation that follows from the covariance of the theory. Alternatively the brane 
may be treated as a common boundary of two distinct regions, Ai^ and A1~ in the 
bulk Ai = 7V1+ UTVl^, which are on the different sides of the brane S = dAi^ = dM.^ . 
The boundary conditions at the brane are given by the Israel equations [50], which 
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correspond precisely to the brane equations of motion. These two approaches are phys- 
ically completely equivalent because the theory is completed with the inclusion of the 
Gibbons-Hawking boundary terms [51], which properly covariantize the bulk Einstein- 
Hilbert action in the presence of a boundary. As a result, varying with respect to the 
metric gives the correct boundary equations as well as the correct bulk. The simpler 
5-function form of the field equations then corresponds to the unitary gauge, realized 
by going to brane Gaussian-normal coordinates, which essentially describe the brane's 
rest frame in the bulk, and then gauge fixing residual gauge invariance. 
The dynamics of the model can therefore be derived from the action 

S = Ml I d^Xy/^R+2Ml f d'^Xy^AK+ f d^Xy/^{Min-a + CmaUer) (2.1) 

Jm Jt, Jt, 

Here gab is the bulk metric with the corresponding Ricci tensor, Rab (in TW = 7\/l+ u 
Ai^). The induced metric on the brane is given by 'jab and its Ricci tensor is TZab, while 
a is the brane tension. The extrinsic curvature of the brane is given by Kab = —^'Cn'jab, 
where Cn is the Lie derivative of the induced metric with respect to unit normal, n", 
oriented from M~ into A4~^; AKab = K^^^ — K~^ is the jump of Kab from ^A~ to Al"*", 
and Cmatter IS the Lagrangian of brane localized matter fields, with vanishing vacuum 
expectation value, because the brane vacuum energy was explicitly extracted as tension. 
In what follows we will use different gauges for the bulk geometry, because the brane 
Gaussian-Normal gauge is very convenient for counting up the modes in the spectrum 
of the theory, while other gauges may be easier to compute the effective actions for 
particular modes. Thus, thinking of the solutions geometrically as a bulk in which the 
brane moves, we will write the field equations which follow from (2.1) as separate bulk 
and brane equations of motion respectively. These are valid in an arbitrary gauge, and 
may be thought of as a breakdown of the full set of field equations on a space with a 
boundary, where the boundary conditions describe a codimension-1 brane. The bulk 
equations of motion are simply the vacuum Einstein equations. 

Gab = Rab — -^Rgab = , (2.2) 

whereas the brane equations of motion are given by the Israel junction conditions [50], 

Qab = Ml A [Kab - K^ab] + MJiUab - ^7^7,b) + '^lab = ^Tab • (2.3) 

^-b - -^^ q::^. ' ^2.4) 

explicitly does not include the brane energy-momentum, 0"7ab. 



where 
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In most of what follows we will impose Z2 orbifold symmetry about the brane 
(section 3.5 will deal with general perturbations). In other words, we will identify A^"*" 
and Ai~ , restricting the dynamics to the Z2 symmetric action given by 

S = 2Ml f (fx^/^R + AMl f d*x^f^K+ f d^ x^/^ {Mill- a + Chatter) (2.5) 

where Kab = K^^^ = —K^j^. The bulk equations of motion (2.2) are of course unchanged, 
while the brane equations of motion simplify to 

Qab = 2M| [Kab - K^ab] + M^iUab ' -7^7a^>) + ^lab = -Tab • (2.6) 

2.1 Background solutions 

Cosmological DGP vacua describe tensional branes in 5D locally Minkowski patches 
glued together such that the jump in extrinsic curvature matches the tension and 
the intrinsic Ricci curvature contributions as in Eq. (2.6). The solutions can be easily 
constructed by taking a bulk geometry which solves the sourceless bulk equations (2.2), 
and slicing it along a trajectory [t{T),R{T)) which solves (2.3). Then R becomes the 
cosmological scale factor and r the comoving time coordinate. Such techniques have 
been used before in the RS2 framework [52, 53]. When the brane only carries nonzero 
tension, its worldvolume is precisely a 4D de Sitter hyperboloid representing the 4D 
de Sitter embedding in a 5D Minkowski space as required by the symmetries of the 
problem [54]. This solution generalizes the geometry of Vilenkin-Ipser-Sikivie inflating 
domain walls in 4D [55], and was in fact also found in [56] in the context of finite 
thickness domain walls. 

In conformal coordinates x" = {x^,y), the full background metric is given by 

ds"^ = Qabdx^'dx'' = a^{y) {dy^ + ^^ydx^dx"") , (2.7) 

where 

^^^dx^'dx" = -df + e^^* dx'^ . (2.8) 

and 

a{y)=exp{eHy), e = ±1. (2.9) 

The bulk spacetime, Ai, is the image of the line < y < oo, with the brane positioned 
aX y = 0. In DGP brane induced gravity theory there exist two distinct branches of 
bulk solutions, labelled by e = ±1. The solution with e = — 1 is commonly referred 
to as the normal branch whereas the solution with e = +1 is referred to as the self- 
accelerating branch, a terminology which will become transparent shortly. The brane 
metric in (2.8) represents the 4D de Sitter geometry in spatially fiat coordinates, which 



cover only one half of the AD de Sitter hyperboloid. The complete cover with global 
coordinates involves the metric ds'^ = —dr"^ + -^ cosh^(ifr)dr23 describing a sequence 
of spatial spheres 5*^, of radius -^ cosh(ift) and spatial line element dfis, which initially 
shrink from infinite radius to radius 1/ H, and then re-expand back to infinity. 



M\nko\Nsk\ 

bulk __ X normal 



de Sitter 
brane 




Figure 1: Embedding of a de Sitter brane in a flat 5D bulk. The brane world volume is 
the hyperboloid in the Minkowski bulk. The normal branch (e = —1) corresponds to keeping 
the interior of the hyperboloid, and its mirror image around the brane. In contrast, for the 
self- accelerating branch (e = +1), we keep the exterior, and its reflection. The latter scenario 
includes the inflating tensionless brane solution. 

The intrinsic curvature H is given by the tension, as dictated by the brane equations 
of motion (a.k.a. brane junction equations) (2.3) at y = 0, 



SM^H^ 
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The solutions are 



H 






l±Wl + 



Mja 
6M| 



(2.10) 



(2.11) 



This equation suggests that there are in fact four possible values of the intrinsic curva- 
ture. However this is not the case. It is easy to see that only two of these solutions are 



-9- 



independent. Indeed, note that a bulk reflection z -^ —z and a time reversal t —^ —t 
map two of the solutions (2.7,2.8) with H < onto the solutions with H > simulta- 
neously reversing the sign of e. Thus without any loss of generality we flx the signs by 
requiring that a positive tension corresponds to positive intrinsic curvature H, so that 






4 



, , Mia 



(2.12) 



We can embed these solutions in the bulk as in flgure 1 [33]. For e = 1, or "self- 
accelerating" branch, we retain the exterior of the hyperboloid. For e = —1, or the 
"normal branch" , we keep the interior of the hyperboloid. It is now clear whence the 
terminology: on the self-accelerating branch, even when the tension vanishes, a = 
0, the geometry describes an accelerating universe, with a non-vanishing curvature 
{H = 2M^/Ml) produced solely by the modiflcation of gravity. The scale of the 
curvature needs to be specially tuned to the present horizon scale of ~ 10~^^ eV, which 
corresponds to about M5 ~ 40 MeV [28, 29], but once this is done one may hope 
to explain the current bout of cosmic acceleration even without any Standard Model 
vacuum energy. The self-accelerating branch of solutions are a distinct new feature of 
the DGP model, they do not exist on Z2-symmetric brane without the induced gravity 
terms on the brane [54, 28]. However, related solutions may arise in theories with 
asymmetric bulk truncations [48, 49]. 

2.2 How do we obtain AD gravity in the DGP model? 

A crucial question is: given the cosmological DGP vacua reviewed above, how could 
there be a low energy 4D gravitational force between masses inhabiting them? Unlike in 
RS2, for solutions given by (2.8) and (2.10), the 'apparent' warping of the bulk cannot 
play a signiflcant role in manufacturing AD gravity at large distances. In RS2 bulk 
gravitational effects pull the KK gravitons away from the brane, strongly suppressing 
their couplings to brane localized matter. As a result, the extra dimension is hidden. 
That does not happen here because the bulk in (2.8) is locally flat. Moreover, on 
the self-accelerating backgrounds the bulk volume is infinite, and so the 4D graviton 
zero mode is decoupled: it is not perturbatively normalizable, and the mass scale 
which governs its coupling diverges. Although the bulk volume for the normal branch 
solutions is finite for finite 1/H, and there is a normalizable graviton mode, its coupling^ 
is Qq oc H/M^, and so it also decouples in the limit H ^ [33]. In fact, from the general 
embedding of a 4D de Sitter hyperboloid in 5D Minkowski space (2.8) we see that the 



''This formula is precisely the analogue of the Gauss law relation between bulk and effective AD 
Newton's constant in models with large extra dimensions [12]. 
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H ^ limit corresponds to taking the radius of extrinsic curvature of the hyperboloid 
on the normal branch to infinity, de facto pushing it to the spatial infinity of Minkowski 
space. In this limit the bulk volume between the brane and the horizon diverges, which 
is why the zero mode graviton decouples. This agrees with the perturbative analysis of 
the H = case of [19, 27] where the zero mode graviton was completely absent. Hence 
4D gravity ought to emerge from the exchange of bulk graviton modes. 

Suppose first that the graviton kinetic terms reside only in the bulk. In an infinite 
bulk, a typical bulk graviton sourced by a mass on the brane will not venture too far 
from the brane because it would cost it too much energy. Nonetheless if kinetic terms 
reside only in the bulk, a typical bulk graviton would still peel away from the brane 
and explore the region of the bulk around the brane out to distances comparable to 
the distance r between the source and a probe on the brane. The momentum transfer 
by each such virtual graviton to the brane probe would be ~ 1/p, where p is the 4D 
momentum along the brane, as dictated by the 5D graviton propagator and brane 
couplings. Thus the gravitational potential would scale as 1/r^, and the resulting 
force as 1/r^. Such force-distance dependence would reveal the presence of the extra 
dimension. This would remain true even when H ^ on the normal branch. Although 
a zero mode is present in this case, it cannot conceal the extra dimension because it 
would still be too weakly coupled to provide the dominant contribution to the long 
range force at sub-horizon scales. 

The induced curvature terms on the brane change this in DGP. In order for this 
trick to work, one needs M4 to be big. In this case, the brane localized kinetic terms 
effectively pull the zero mode gravitons closer to the brane, making their exploration 
of the bulk at distances shorter than Vc ~ M'l/2Ml energetically costly [19, 27]. This 
alters the scaling of the momentum transfer to 1/p^ for momenta p > M^/M^, which 
in turn produces a force which scales as 1/r^. This is manifest from the explicit form 
of the graviton propagator projected on to the Minkowski brane (i.e. the H = limit 
of the normal branch solutions of (2.8) and (2.10)) [19, 27]: 

From the AD point of view, the graviton resonance which is exchanged is composed of 
massive tensor modes, and so it will contain admixtures of longitudinal gravitons. This 
is encoded in the propagator (2.13) in the coefficient 1/3 of the last term of the spin 
projector, as opposed to 1/2 which appears in the linearized limit of standard 4D GR. 
This difference is an example of the venerated Iwasaki-van Dam-Veltman-Zakharov 
(IvDVZ) discontinuity of modified gravity [57], and signifies the persistence of a scalar 
component of gravity in the theory, that could confiict with the known tests of GR. 
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However, it has been argued for massive gravity [58] and similarly for the DGP model 
[59] that the extra scalar may be tamed by nonlinearities once the correct background 
field of the source is included. The idea is that the perturbative treatment of the scalar 
graviton breaks down at a distance scale ry first elucidated by Vainshtein [58]. For 
DGP, for a source of mass m, this new scale is given by ry ~ {mr'^/MjY^^ [60, 61]. 
Below that distance, one can't trust linearized perturbation theory and must re-sum the 
background corrections, which should presumably decouple the scalar graviton mode. 
Similar weakening of the scalar graviton coupling may occur at cosmological scales if 
the universe is curved. 

This scale dependence of the scalar graviton couplings has very interesting and im- 
portant implications for the DGP setup. It has been pointed out [35, 36, 37, 38] that the 
effective field theory description of DGP gravity will suffer from a loss of predictivity 
due to the problems with strong couplings at distances rstrong ~ {r^/M^y^^, which could 
be much larger than the naive UV cutoff. The most recent analysis of this issue [38] 
however suggests that the brane nonlinearities may push the scale of strong coupling 
down, to about fstmng ~ ^strong / \/ Mearth / M4 ~ Icm ou the surface of the Earth, pos- 
sibly making DGP marginally consistent with current table top experimental bounds 
[6, 7]. In what follows we will assume this claim [38] and imagine that we work in the 
perturbative regime of DGP, although we feel that this issue deserves further atten- 
tion. We note that the exploration of DGP with gravitational shock waves [33] shows 
that the scalar graviton decouples from the background of relativistic sources, indicat- 
ing that the coupling is effectively suppressed by the ratio of \/iT^\iy/T^^\iu of the 
source. Note, that this is not enough to ascertain that a theory is phenomenologically 
safe. For example, a Brans-Dicke theory will admit identical shock waves as ordinary 
GR for any value of the Brans-Dicke parameter u, while the observations require that 
u > 5000. Thus one still needs to study the model for slowly moving sources to check 
if the predictions agree with observations. However one may hope that the strong cou- 
pling problems might be resolved in a satisfactory fashion. After all, the shock waves 
[33] remain valid down to arbitrarily short distances from the source, behaving much 
better than they are entitled to given the concern about the strong coupling problems. 

In what follows we will focus on uncovering the ghosts (and/or other instabilities) 
on the self-accelerating branch, and a discussion of their implications for DGP. Before 
we turn to this, we should stress that there is no technical inconsistency between our 
results and the earlier claims that there are ghost-free regimes in DGP [19, 27]. Indeed: 
starting with the backgrounds of the family (2.8,2.10) and fixing M5 and M4, the only 
way to consistently take the limit to // — ^ is to pick the normal branch solutions 
and dial the brane tension a to zero. In this way one reproduces the H = brane 
backgrounds with fixed M5, M4 that were studied in [19, 27]. Moreover, ghosts may 
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also be absent if the brane geometry is anti de Sitter, as opposed to dS [62]. Thus 
the results of the perturbative analysis of [19, 27], implying the absence of ghosts and 
other instabilities on H = branes, applies only to the normal branch backgrounds 
of DGP (2.8), (2.10). In fact, our results will confirm this for the general H ^ 
backgrounds of the normal branch, showing that they are ghost-free. However the 
analysis of [19, 27] has nothing to say about the self-accelerating branch solutions, and 
specifically about the o" = limit, that describes a universe where cosmic acceleration 
arises from modification of gravity alone. In what follows we will confirm that in all 
those cases there are ghosts, which invalidate the self-accelerating branch solutions in 
their present form as realistic cosmological vacua. 

3. The occult sector of DGP 

We now turn to the exploration of the spectrum of light modes in the gravitational 
sector of DGP, around the cosmological vacua (2.8), (2.10). We will confirm that there 
are ghosts in the 4D effective field theory description on the self-accelerating branch 
of DGP solutions. More specifically: in the 4D effective field theory which describes 
the perturbative regime of self-accelerating branch of DGP backgrounds (2.8), (2.10) 
between the Vainshtein scale ry and the scale of modification of gravity r^ = 'mij2Ml 
there are scalar fields with negative, or vanishing, kinetic term around the vacuum, 
which couple to the brane-localized matter with at least gravitational strength. Now, 
this may appear surprising at the first glance: there are no ghosts in the action (2.1) 
of the full 5D bulk theory. Indeed, the full bulk Lagrangian in (2.1) does not appear 
to contain any instabilities. However, the background solutions (2.8), (2.10) of (2.1) 
involve an end-of-the-world brane, which is a dynamical object, whose world-volume is 
determined by (2.3). The problems arise because the brane will curl up and wiggle when 
burdened with a localized mass, in a way that alters the gravitational fields of the source 
mass, spoiling the AD guise of the theory. Thus the perturbative ghost encountered 
in AD theory is really a diagnostic of the failure of the AD perturbation theory to 
describe the dynamics of the long range gravity on the self-accelerating solutions. Thus 
although the applications [28, 29, 63] of the self-accelerating solutions to AD cosmology 
are interesting and tempting, the presence of the ghost renders them unreliable at the 
present stage of understanding of the theory, and hence de facto inadequate as a method 
of accommodating the present stage of cosmic acceleration. 

In the following subsections we will identify the independent degrees of freedom 
describing small perturbations around DGP vacua in both branches, derive their lin- 
earized equations of motion and solve them. We will then compute the four dimensional 
effective action, isolate the ghost of the AD theory, and discuss its consequences. 
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The physical interpretation of these solutions is based on the mathematical analysis 
of a differential operator derived by considering perturbations of Einstein's equations: 
the Lichnerowicz operator Aihab- This operator acts on a five-dimensional spacetime 
with a timelike boundary (the brane). We can solve these perturbation equations in 
whatever gauge we like, however, in order to get a braneworld interpretation of the 
results, the cleanest procedure we can follow is to separate this problem (operator 
plus space on which it acts) into a direct sum of a purely four- dimensional operator 
acting on a four- dimensional spacetime, and a self-adjoint ordinary differential operator 
acting on the semi-infinite real line. Obviously this latter operator acts on the space 
perpendicular to the brane, and hence to really benefit from this factorization, in these 
coordinates the brane should be held at a fixed coordinate position. Once we have 
made this decomposition, we will be able to identify the physical states and their norms 
from the braneworld point of view. To this end, we should write the perturbation in its 
irreducible components with respect to the braneworld, correctly identify the degrees of 
freedom corresponding to "motion" of the brane, and reduce our perturbation equations 
to a self-adjoint form. 

3.1 Learning to count: mode expansion 

We turn to the linearized perturbations hab{x,y) about the background metric (2.7), 
(2.8), (2.10), defined by the general formula 

ds^ = a\y)(^%b + a{y)-^^^hab{x,y)yx''dx\ (3.1) 

where we use the shorthand ^ahdx"'dx'^ = dy"^ + '^^^dx^dx" . Note that aiy) = exp{eHy) 
as specified in (2.7), (2.8), (2.10). From now on, we will raise and lower 4D indices 
(/i, i^, ■ ■ ■) with respect to the AD de Sitter metric •y^^, and designate 4_D de Sitter 
covariant derivatives by D^. Our normalization convention for the perturbations in 
(3.1) refiects after-the-fact wisdom, in that it simplifies the bulk mode equations to a 
Schrodinger form, as we will see later on. 

Since the spacetime ends on the brane, if we fix the gauge in the unperturbed 
solution (2.7), (2.8), (2.10) such that the brane resides at y = 0, a general perturbation 
of the system will also allow the brane itself to fiutter, moving to 

y = F{xn . (3.2) 

The explicit expressions for the perturbations hab, F are obviously gauge-dependent. 
Now, to consider their transformation properties under diffeomorphisms 

y -^ y' = y + C{x,y), 

x^ -^ x'^ = x^ + x^{x, y) , (3.3) 
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we should first classify them according to different representations of the AD diffeonior- 
phism group as 

h^i, , a tower of 4D tensors ; 

J hya , a tower of 4D vectors ; 
perturbations = < , r .t^ i (3-4) 

' hyy , a tower or AD scalars ; 

F , a single 4D scalar. 

comprising in total 10 tensor + 4 vector + 1 scalar towers = 15 towers of degrees of 
freedom plus one more AD scalar, i.e. precisely the number of independent fluctuations 
of a symmetric 5x5 bulk metric and the brane location. Clearly, not all of these degrees 
of freedom are physical: some can be undone by diffeomorphisms (3.3). Indeed, we can 
easily derive the explicit infinitesimal transformation rules. 



h'yy = hyy-2a'/'(^dyC + eHc), 
F' = F + Cy=o, (3.5) 



where we have used that dya = eH a. 



In order to have a clear braneworld interpretation of variables, we find it convenient 
to work in the Gaussian-normal (GN) gauge (see e.g. [64]), in which any orthogonal 
component of the metric perturbation vanishes. Given any perturbation (3.4), we can 
transform to a GN gauge by picking the gauge parameters (, x^ 



"yy •> 



Xm = / dya'^^^hy^- Df, dyC, (3.6) 

which set h' and h' to zero. This still leaves us with 10 components of /i^,y and the 
brane location F (omitting the primes), accompanied by 5 residual gauge transforma- 
tions 



X" = xS(a^) + -^D-f(x) , (37) 

tha 

which can remove several more mass multiplets from the perturbations. However these 
could only be zero modes of some of the bulk fields, because of the restricted nature 
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of the bulk variation of (3.7). Rather than completely gauge fix the perturbations 
now, it is more useful to resort to dynamics to find out which of the modes h^^^ F 
are propagating and which are merely Lagrange multipliers. To this end we can first 
decompose the tensor h^u in terms of irreducible representations of the diffeomorphism 
group. This yields (for a proof, see Appendix (6)) 



V = h^Hu + D^,A^ + D^A^ + D^D^(I) - -^^^D^cj) + -7^.^ , (3.^ 



where /i^]^ is a transverse traceless tensor, D^hT^ ^^y = hP- ^^ = 0, with 5 components, 
A^ is a Lorentz-gauge vector, D^A'^ = 0, with 3 components, and and h = h'^^ are 
two scalar fields (such that they correctly add up to the total of 10 degrees of freedom). 
To get some feel of the dynamics before looking directly at the field equations, 
we can check how these modes transform under the residual gauge transformations 
(3.7). Substituting the residual gauge transformation (3.7) into (3.5), we find that 
the surviving, symmetric, tensor mode in the GN gauge and the brane location field 
transform as 

F' = F + f, (3.9) 

where Ofj,u = D^Du + H'^^^,^, and we have used that a(0) = 1. If we further split up 
the gauge transformation parameter Xoa* = ^^ + I-^m^' where D^£^ = 0, and apply 
the decomposition (3.4) of h^^, into the irreducible representations /i" , A^, and h to 
(3.9), after a straightforward computation we find that the irreducible representations 
transform according to 

r /TT _ r TT 
"'fJ.U "'IJ.U ! 

eH 

h' = h- a^'^D^uo - '^D^f - 8eHa'/^f , 

eH 

F' = F + f. (3.10) 

Note that while the decomposition (3.8) of a general hf^i, into irreducible represen- 
tations of the diffeomorphism group is kinematically unique, implying the breakdown 
of the residual gauge transformations as in (3.10), it does not - in general - guarantee 
that different modes won't mix dynamically. Indeed, in writing (3.8) we are implicitly 
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assuming that different irreducible transformations live on different mass shells, and 
hence cannot mix dynamically at the quadratic level. This can be glimpsed at, for 
example, by noting that while the symmetries of the problem allow us to write the 
couplings like hJl^,{giD'^D'^ + g2l^^)(t> etc, the TT conditions for /i" would imply that 
such couplings are pure boundary terms for non-singular couplings gi,g2- While this is 
true in general, the situation is considerably subtler when the representations become 
degenerate. In this instance the decomposition (3.8) requires more care. New acciden- 
tal symmetries mixing different representations, notably tensor and scalar, may arise, 
modifying (3.10) and dynamically mixing the modes. This occurs in the vanishing 
brane tension limit on the self-accelerating branch of DGP. We will revisit this limit in 
more detail later on. 

Keeping with the general situation for now, the transverse-traceless tensor hj^^, is 
gauge invariant, while the vector and the scalars are gauge dependent - we can gauge 
away the zero mode of the vector and one of the scalars. In addition, we see how 
the motion of the brane can be gauged away, choosing / = —F to set the location of 
the brane to y = 0. By doing this, we are explicitly choosing coordinates which are 
brane-based, and the metric perturbation (3.8) has an explicit O^^F term describing 
the brane fluctuation. In the brane-based approach, we have completely and rigorously 
separated the Lichnerowicz operator into brane parallel and transverse parts. However: 
once we have taken these coordinates, we do not have the liberty of making residual 
gauge transformations parameterized by / in (3.10), because they would move the 
brane from the coordinate origin. In effect, the brane is tied to the dynamical fields 
and h in the bulk, but its fluctuation F turns into a Goldstone boson of the system. 
We emphasize that this is a gauge choice. We are choosing the brane-GN gauge to 
make the separation of the Lichnerowicz operator mathematically clean. However, one 
can also choose to allow the brane to fluctuate freely (and indeed the effective action 
computation is better done this way) by having a bulk-GN gauge, with the /-gauge 
freedom in (3.10), and the brane sitting aX y = F. In this case, there are no flxed 
terms in the perturbation, and the brane motion enters into the boundary condition 
via evaluation of the background solution at y = F. The actual equations of motion 
and boundary terms in both gauges are identical, giving the same physical results and 
the same dynamical scalar flelds. Thus, explicitly in brane-GN gauge: 

V = C + ^M^- + ^-^M + (^M^ - i^^V) <P + -^^>^-^ + i^V • (3-11) 

To proceed with setting up the problem, we derive the fleld equations for the ir- 
reducible modes. Having restricted to the family of brane-GN gauge perturbations 
(3.11), we can substitute them in the fleld equations (2.2), (2.3) and after straightfor- 
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ward algebra write the linearized field equations in the bulk, 

5Gab = , 



where 



a^^^SG^, = X^,{h) 






9i£ 

4 



( V - ^It^y) 



a"HG 



fj-y 



d_ 
dy 



3eH 



D" {Ku - /i7, 



^JLU ) 1 



t'lHG, 



3eH 



yy 



d_ 
dy 



3eH 



h-- ID^D' 



flUj 5 



7^^(1)2 + 3H^)] h 



tJ.U ) 



(3.12) 

(3.13) 
(3.14) 

(3.15) 



and 



7^, [D^D^K^ - {D^ + H^) h] 



1- 

2 



and on the brane, 



'/iZ^ 



d_ 
dy 



3eH 



{h^,v - hj^i,) 



It 

2 1^" 



(3.16) 



(3.17) 



y=0 



Before we proceed with the details of the mode decomposition of this system by direct 
substitution of (3.8), we note that the Lorentz-gauge vector A^ turns out to be a free 
field in the linearized theory in flat bulk. Thus the solutions for A^ decouple from 
the brane-localized matter in the leading order. They are irrelevant for the stability 
analysis which is our purpose here. Hence we will set A^ = from now on, assuming we 
have arranged for bulk boundary conditions which guarantee this in the linear order. 

3.2 Fluctuations around the vacuum 

First note that independently of matter on the brane, the yy and yn equations must 
be identically satisfied. In conjunction with the trace of the /iz/ equation, this can be 



seen to imply that a gauge can be chosen in which h = 0, and O^ 
we have no matter on the brane, then we see that 



0. If in addition 



O^xF = (D^ + AH^)F 



0, 



(3.18) 



and so the metric perturbation (3.11) is completely transverse and traceless. 
The remaining ^v bulk equations then simplify considerably to give 



D' - 2H^ + 



dy^ 



h^,y{x,y) = 0. 



(3.19) 
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with the boundary condition 

= (3.20) 



M! {D' - 2H') V + 2Ml ^-^ - ^) V 



y=0 

Now, to solve this equation we should carefully decompose the tensor into orthog- 
onal modes, which in general do not mix at the linearized level. Those are exactly the 
irreducible representations we discussed previously. Thus using linear superposition, we 
can expand the general metric fluctuation in /i" and 0, the latter of which couples to 
the field F through the boundary condition (3.20), leaving the TT-tensors with entirely 
homogeneous boundary conditions. We therefore write 

V = E^™(^)4i^^(^) + ^it^ + ^^^-^ (3.21) 

m 

where we have performed the mode expansion W-^^{x,y) = 'Yl,m'^-m{y)Xjj^ {x) ■, in terms 
of the AD modes Xm™ (^) which satisfy {D'^ — 2H'^)x^lu = rr?X]Z . We have also 
defined the scalar mode /i^tj/ = 0^y(\)^ and separated variables in the scalar field by 
setting = W{y)(^[x\ where is a general AD tachyonic field obeying 

(L)^ + 4ij2)0 = . (3.22) 

This tachyonic mode is present whenever we compactify the theory on an interval 
with de Sitter boundary branes. It can be traced back to the repulsive nature of 
inflating domain walls [55]. Here, it is simply an indication that a multi-de Sitter 
brane configuration requires a special stabilizing potential, as is familiar already in the 
context of RS2 braneworld models [65, 66, 67]. This kind of an instability is generically 
much slower and hence less dangerous than the ghost, as it is governed by the scale r ~ 
1/mtachyon ~ H ^ that is as long as the age of the universe. When tension is positive, 
this mode therefore remains largely harmless for the phenomenological applications of 
the theory. However, on the self-accelerating branch it mixes with the ghost in the 
vanishing tension limit on the self-accelerating branch, and becomes the ghost itself for 
negative tension, as we will see later. 

We now turn to the analysis of the TT perturbations which form the main part of the 
propagator, and determine the norm on the transverse y-space. The bulk field equation 
(3.19) and the boundary condition (3.20) reduce to the boundary value problem 

^'m{v) + I "^^ T- ) Urn{y) = , 



Mi 



u'miO) - ^Mm(O) 



+ |m2M|M„(0) = , (3.23) 
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which is self-adjoint with respect to the inner product 



{u\v) 



dy {Ml+Mi6{y))u{y)v{y) 



2M| 



dyu{y)v{y)+Miu{0)v{0). (3.24) 



The eigenmodes Um with different eigenvalues m are orthogonal. We choose the normal- 
ization such that the discrete modes, if any, satisfy («„!«„) = Smn, while the continuum 
modes satisfy {um\un) = 5{m — n). This is simply a reflection of the fact that far from 
the brane the bulk modes behave just like bulk plane waves, and the AD mass is pre- 
cisely the Pj^-component of the bD momentum. 

To determine the spectrum of the boundary value problem (3.23), (3.24) we rewrite 
the boundary value problem (3.23) as a Schrodinger equation 

^4- + A,n^ - 3mM 



" -L 



ra 



■Ml 



Ur 



0. 



(3.25) 



It is now clear that the solutions of (3.23) must fall into two categories: {%) one dis- 
crete mode for each branch, localized to the (5-function potential, if it is normalizable 
according to (3.24), and {ii) a continuum of 'free' modes, gapped by m > |i/. 

• m^ < ^^-: the normalizable solution of (3.23) in the bulk, representing a single. 



light, localized graviton on each branch, with a mass 



niA 



Ml 

Ml 



3H 



2M| 



:i+e) 



(3.26) 



fixed by the boundary conditions (3.23), and wave function 

"3M|if-2M|(l 



Umiy) = amexp{-\my) 



a. 



1 

M4 



3M|iJ - 2M|e 



where \r 



4 



— m? and {u\u) = 1. On the normal branch (e 
mrf = 0. 
On the self accelerating branch (e = +1), 

< m^ < 2i/2 ^ for (7 > ; 



(3.27) 



(3.28) 



ml > 2H^ 



for o- < . 



(3.29) 



Herein is our first glimpse of the tensor ghost: for positive tension, the localized 
light graviton mode on the self accelerating branch lies in the forbidden mass 
range < m'^ < 2H^ discussed in [44, 45, 46]. Its helicity-0 component is the 
ghost, as we will review later on (see Appendix (7)). 
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^2 y 9H.^ ^YiQ (5-function normalizable modes are 



772 
Um{y) = ttm Sm{uJmy + (^m) , Cim = \ / ^TTg , (3.30) 



where Um = J rn? — ^^ and {um\urh) = ^(^ — f^)- The integration constant 5m 
which solves the boundary condition (3.23) is 

Turning now to the scalar component hifj{x,y) = W{y)0^,y(f){x), it is not difficult 
to see that it obeys 

{D' - 2H') hl^J = 2H%^J , (3.32) 

(equivalent to a AD mass m^ = 2H^). The bulk equation (3.19) then yields the wave 
equation for W, 

W"iy)-^Wiy) = 0. (3.33) 

The boundary condition (3.20) enforces a relation between and F: 

(,r,0,-(5..-^),.,0,),^.(.-.^).. ,3.34, 

The wave function solutions for either of the DGP branches are 

Wiy) = aexp (-^v) + /^exp (^y\ . (3.35) 

From (3.24), the norm is determined by L dyW'^{y), where the lower limit of inte- 
gration accounts for the unperturbed location of the brane at y = 0, around which we 
impose the Z2 symmetry. Thus the a-mode is normalizable but the /9-mode is not. 
We therefore set f3 = 0. This choice, at least in principle, corresponds to prescribing 
boundary conditions at infinity, which ensure the brane is an isolated system. Thus 
setting (3 = 0, and separating the variables by setting (p = F in (3.34) we find 

« = 2 1-— tt:^ • (3.36) 



Ml J V Ml 



However, we must be mindful of this choice because of the possible interplay with the 
brane bending term (3.10), as we will see next. 

Now: on the normal branch (e = — 1), it follows from (3.10) that the normalizable 
a-mode is gauge-dependent: in fact, it is of the same form as the brane-bending mode 
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since it is proportional to a^/^ = exp (— -f !/)• On the other hand the non-normahzable 
/9-niode is gauge-invariant by itself, and so setting it equal to zero is straightforward. 
Then (3.36) gives a = 2/H, which means that the brane boundary condition (3.20) in 
fact precisely sets the normalizable gauge-invariant mode a(f) — 2F/ H to zero. Hence 

Ku = hj^^ (3.37) 

Thus the net effect of the a-mode is to undo the brane bending. This is because 
the translational invariance of the brane-bulk system, which yields the residual gauge 
symmetry (3.10) is linearly realized in the presence of the brane, which imposes gauge- 
invariant boundary condition, so that the normalizable bulk mode and the brane bend- 
ing completely compensate each other. Put another way, the only consistent matter-free 
solution for the normal branch DGP brane is where the brane does not move from y = 0, 
and only TT GN perturbations in the metric are allowed. This, of course, should have 
been expected all along, as it is just the statement that the radion field decouples in 
the case of single UV brane with 4D Minkowski or de Sitter geometry embedded in the 
standard way in 5D Minkowski or AdS space. Here we see explicitly how gauge invari- 
ance and normalizability enter this subtle conspiracy to remove this mode, essentially 
allowing that any scalar bulk perturbation localized to the brane can be bent away. 

On the self-accelerating branch (e = +1), the situation is very different: now, the 
normalizable scalar mode is gauge-invariant by itself. The non-normalizable /3-mode 
is not, and so imposing boundary conditions which require (3 = breaks the residual 
gauge invariance (3.10). The brane bending mode F is the Goldstone field of the broken 
symmetry, and the brane boundary condition (3.20) for a generic value of H (i.e. for 
non-zero tension) yields 

2 



" = -^ 



Ml - MjH 
2M| - MjH 



(3.38) 



which pins the Goldstone F to the normalizable gauge-invariant scalar perturbation 0: 

e-^y'^O^^F. (3.39) 



/i(<^) = - — 



Ml - MlH 



2M| - MjH 



This perturbation represents a genuine radion, or physical motion of the brane with 
respect to infinity. Although our choice of brane-GN gauge fixes the brane to the 
coordinate position y = 0, it does so at the cost of, this time, breaking the residual 
gauge symmetry generated by / in (3.10) and introducing the explicit "book-keeping" 
On„F term in /i^j,, which is the remnant of the translational zero mode of the brane. 
Had we instead allowed the brane position to be arbitrary, aXy = F, (without the O^i^F 
term in (3.11)), the boundary conditions aX y = F would still have had the same form. 
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since the F-terms would have entered when evaluating the background at nonzero y. 
Both approaches are completely equivalent, the former being more suitable to a brane 
based observer and the latter to an asymptotic observer. The gauge transformation 
between these is a y-translation, which therefore corresponds to real motion of the 
brane, just as in the 2-brane RS case [41]. The absence of this mode on the normal 
branch reflects the fact that there is no distinguishable motion of an individual Z2 
symmetric brane. 

When the tension is different from zero, the solutions Xm™ are precisely the TT- 
tensors /i" of (3.4) from the previous subsection. The scalar mode hilJ has eigenvalue 



(m) 



m = 2H , as seen from (3.32), and the eigenvalues of the eigenmodes xl^v are all 
different from 2H^ when a 7^ 0. Thus the scalar mode 0, disguised as the tensor h 



.M w^v.^^ ..M 






is orthogonal to all Xm^^ • Hence Xm^^ coincide with the TT tensors /i", and so when 



/ii^' 



there is no matter on the brane, the solutions are given by 

/■oo 

V(x,i/) = a,n,e-^™<^^x!.7Ha^) + /^ dm M„(i/)x!.™Ha:) 

This solution clearly remains valid on the normal branch even in the limit of vanishing 
tension, a — * 0, and for the full range of o" < 0, because when e = — 1 the potentially 
dangerous O^^F terms vanish identically. 

However on the self-accelerating branch where e = +1 the solution (3.40) - as it 
stands - fails when the tension vanishes, o" = 0, because of the pole in 0, or a, (3.38), 
(3.39). Indeed, (2.12) implies that when a —^ 0, H ^ 2Ml/Ml, and so the parameter 
a in (3.38) diverges. Thus the mode as given by (3.39) is ill-defined in this limit. At 
a glance, noting that the coefficient of in (3.34) vanishes, one may interpret equation 
(3.38) as implying F = 0, thus fixing the brane rigidly at y = 0, and allowing to 
fiuctuate independently of F. However, in light or the residual gauge transformations 
(3.10) and our gauge-fixing f3 = 0, that removed the non-normalizable gauge-dependent 
bulk scalar, setting F = also would completely break the residual gauge symmetry 
group. This is dangerous, since it may miss physical degrees of freedom, which warns us 
against such a quick conclusion. To see what is really going on we must tread carefully. 

What's going on when the tension vanishes is that the mass of the localized tensor 
mode on the self-accelerating branch approaches m^ = 2H^, as is clear from (3.26). 
Further, the bulk wave function of the lightest localized tensor (3.27) converges to 
exp(— Am^y) = exp{—Hy/2), i.e. it becomes identical to the bulk wave function of 
the gauge-invariant scalar mode h/Tj . Thus the lightest tensor, h^v ■, and the scalar 
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hjfj = O^yCJ) become dynamically degenerate^ and can mix^ together: they both solve 
the AD field equations [D"^ — 4:H'^)Xfiu = and have formally the same tensor structure. 
Now, it has been noted by Deser and Nepomechie [43] that in the special case when 
the mass of the massive AD Pauli-Fierz theory in de Sitter space equals 2H^, the 
theory develops an accidental symmetry[43, 44, 45, 46]. The tensor dynamics becomes 
invariant under the transformation x\iv -^ x\i.v + O^^^^ where ^ is any solution 
of the equation [D^ + AH'^)tD = 0, affecting only the helicity-0 component of xliiy ■ 

Lifting this symmetry to the present case is considerably more intricate because of 
the degenerate scalar h/jj . Noting first that the wave profile of the lightest localized 
tensor is now e~^^/^ = l/^/a^ the accidental symmetry of [43] shifts the bulk TT-tensor 

by 

h';r.=hll + a-'/'0,.^. (3.41) 

However given the higher- dimensional origin of the perturbations /i^jy we cannot ar- 
bitrarily shift these modes around. The only gauge generators available to us, that 
could in principle generate such shifts, are the residual gauge transformation rules of 
(3.9). However as is clear from (3.9), none of the residual gauge transformations have 
the correct bulk wave profile to yield (3.41). Thus the transformation (3.41) must be 
understood as the Stiickelberg symmetry of the problem: shifting xlii' by a -(9 piece 
must be compensated by shifting another field in the decomposition (3.8) to keep the 
total metric perturbation hfj^i, invariant. The only available mode with the correct wave 
profile, and the correct tensor structure so as not to break the diffeomorphism invari- 
ance, is the normalizable scalar that is invariant under (3.10). Thus the scalar must 
now be promoted into a Stiickelberg field for i}. Note that this is completely analogous 
to rewriting the massive U{1) gauge theory in the Stiickelberg form, formally giving up 
on the Lorentz gauge for A^ by introducing the Stiickelberg scalar field. 

So to properly account for the accidental symmetry on the self-accelerating brane 
in the vanishing tension limit generated by i}, we must enhance the residual gauge 
transformation group (3.10) by also including in it 

C = Kl + a''/^0^,^ , 0' = - a-'/^^ . (3.42) 

Hence once we insert a = explicitly in the field equations (2.2), (3.20), we can 
separate the field equations for the scalar and the lightest tensor from each other only 
after explicitly gauge-fixing the Stiickelberg symmetry generated by i). The full field 
equations are merely covariant under it because of the scalar field 0. Once we have 



^This mixing has been noticed as the resonance instabiHty in the shock wave analysis of [33], and 
discussed at length in [68]. 
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gauge-fixed the brane at y = 0, the boundary conditions at the brane will really relate 
F to a linear combination of the helicity-0 tensor and the gauge-invariant normalizable 
scalar 0. A simple way to think about the boundary conditions is to fix the i) gauge 
by completely removing the helicity-0 mode from the tensor and absorbing it into 0. 
Then the brane boundary condition (3.38) just states that this -(^-gauge fixed field is 
fluctuating freely - but it does not disappear from the spectrum. Indeed, we can go 
to a different -i^-gauge, fixing it now such that the is completely eaten by the tensor 
Xnu , which regains its helicity-0 component. This of course is completely equivalent 
to the unitary gauge of a theory with the Stiickelberg fields, where the gauge fields eat 
Stiickelberg and gain mass. This is crucial for the failure of the ghost to decouple in 
the tensionless brane limit, and is a simple way to understand the analysis of [68]. 

However, neither of these gauges is convenient for the computation of the effective 
action, to be pursued later on. Instead, we can pick another i) gauge by taking the 
general solution (3.40) and defining a i} which removes the pole in (3.40) as o" ^ and 
produces a smooth limit [39, 68]. We can do this by a shift 

c^m.xi^'Hx) = n^uix) - a 0^,F . (3.43) 

Substituting this in (3.40) yields 

poo 

h^u{x,y) = e-^-'^y-H^^uix) + / dm Umiy)x^;:\x) 

Then carefully taking the limit a ^ (noting that a oc 1/a, and Xm^ = H/2 + 0(o")) 
yields 

h,,{x,y) = e-Ty(n,.{x)-yO,,FJ+ / dm uUy)x'p\x) + -e^'y/^O^^F . (3.45) 

where the 4D tensor H^u satisfies ^ 

{D^ - AH^)n^, = -HO^,F . (3.46) 

Note that the trace of this equation yields O'^xF = {D^ + 4:H'^)F = 0. In this gauge, 
a way to think about Eq. (3.46) is to view the field F as the independent degree 
of freedom, and the helicity-0 component of the graviton TY^j, as being completely 



®This follows from the cr ^ limit of (D^ — m^ — 2H^)T-ifj,i, from (3.43), or can be readily derived 
from the equations of motion (3.19, 3.20). 
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determined by the source F. Yet, the brane locahzed matter can only feel its influence 
through the couplings to the helicity-0 component of Ti^y, as can be seen from Eq. 
(3.45), which shows that on the brane at y = the oc 0^i,F terms vanish. In this way, 
the tensor H 



fit/ 



[x) retains five physical, gauge-invariant degrees of freedom precisely 
because of this mixing with F, inherited from the bulk scalar 0. In effect what happened 
in the limiting procedure is that the pole of (3.40) was a pure gauge term of the 
Stiickelberg gauge symmetry, and was absorbed away by the choice of i), leaving in its 
wake the smooth function (3.45). 

3.3 Linearized fields of matter lumps 

Here we include the contributions from localized stress-energy lumps on the brane. In 
the linearized theory, the general solution is a linear combination of the homogeneous 
solution, given by (3.40) or (3.45), describing propagating graviton modes, and a par- 
ticular solution comprising of a TT piece and the relevant brane bending term which 
describe the response of the fields to the source. Thus we write: 






eH 



(3.47) 



where the fields x^u{x) and f{x) are the sought-after particular solutions that include 
the effects of the brane sources. They must be the solution of the boundary value 
problem, 

92 9/J2- 



D' - 2H^ + 



dy2 4 



Xt^u{x,y) =0, 



A. 

dy 



3eH 

2 



XfiU 



y=o 



-2 {Ml - MleH) {O,, - 0\%,) f = -\T^ 

Tracing the boundary condition immediately gives 

T 



flU 



.X 



0\f= [D^ + AH^]f{ 



X] = — 



_12{MI - MjeH) _ 



(3.48) 



(3.49) 



This fixes f{x) completely, as any homogeneous brane-bending term is accounted for 
in h 



{horn) 



liu 



We can write the particular solution as a spectral expansion, using the properties 
of the 4D mass eigenmodes determined in section 3.2. Formally, we consider the 4D 
differential operator D^—2H^ whose tensor spectrum is given by TT-tensors x/S obeying 
{D'^ — 2H'^)x^u = P^XiSv, and expand the solutions and the sources as 



XtJiAx^y) = / Vp{y)x}^Ax) 



-TpX^ 






X] 



(3.50) 
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where 



'^flU\-^ ) 



T^, - 4(M| - MleH) [D^D, - %,{D^ + 3H^)] f{x) 



(3.51) 



The Xf^i^ tensors are an ortho normal basis of the spectrum of D^ — 2H^, whose eigen- 
values j9^ 7^ m'^ are taken to be off mass shell as is usual in the inhomogeneous problem. 
Here, J is a generalized sum, accounting for the integration over the continuum part 
of the spectrum and the summation over the discrete, localized modes. Then the 
boundary value problem (3.48) reduces to the system 



Vpiy)+ (p^ J- ]vpiy) = 0. 



Ml 



<(0) - -^vp{0) 



p'MivpiO) = -ir, 



2 P ' 



(3.52) 



extending (3.23) of section 3.2 with a source term Tp. We can write the solutions 
in terms of the on-shell eigenfunctions Um{y)- First, rewrite (3.52) as 

9H^ ,M^ , „ 1 r, 



VpS 



< + 



P'- 



M^ 
4 +ij^f'-^^H)5{y) 



Ml 



yiy) ■ 



(3.53) 



Then expanding as Vp{y) = Vpm^Um^^y) + J^,^ dmVpmUm{y), substituting in (3.53) and 
comparing with (3.25), and finally using the orthonormality of the eigenmodes Um with 
respect to the inner product (3.24), we find 



My) 



Mmrf(y)Mmrf(0) 

p2 — m^ 



dm 



3H_ 

2 



Um{y)Um{0) 
p2 _ jj^2 



Tp. 



Hence we can rewrite the linearized field of the matter on the brane as 

XtJiu{.x,y) = / rfVv^ G^u"'^{x,y;x',0)Ta(3{x') , 
where the Green's function is given by the eigenmode expansion 

UmSy)UmM 



(3.54) 



(3.55) 



G^,"^{x,y;x',0) 



p^ — rn?^ 



°^ , Mm(l/)Mm(0) 



3-g 

2 



p2 _ ^ 



>)(^)^*(p)«/5(^')^ 



(3.56) 

where the asterisk denotes complex conjugation. Thus in the case when the only source 
of perturbation of the vacuum is matter on the brane, and the propagating geometric 
modes are not excited, the brane geometry is perturbed by 



2e 

h^y = Xf^u{x, 0) + —o^^f 



(3.57) 
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where / is given by (3.49) and x^lu{x,0) by (3.55). Note however that one must treat 
the Green's function (3.56) with care, because the summation ^ over the continuum 
has a branch cut at m'^ = 9iJ^/4, which can be seen from the form of the continuum 
eigenfunctions presented in section 3.2. 

The ghost is hidden in the locahzed mode contribution to G^^°''^{x,y;x',0), (i.e. 
the "'"d y "™d — term), specifically, it resides in the helicity-0 component. We could 

p m^ 

divine the ghost by computing the residues at the pole p^ = m^ of the propagator. 
Alternatively, as we will do in the next section, we can simply compute the effective 
action for small metric fluctuations and unveil the ghost-like behavior from the negative 
contributions to it. 

3.4 Forking the ghost: calculating the effective action 

Let us now fork^ the ghost: we compute the effective 4D action of normalizable small 
perturbations considered in the previous section, that will serve as a straightforward 
diagnostic of the ghost. We start with the general case of non- vanishing tension, a 7^ 0, 
and consider the limit a = separately. Let us consider the general metric perturbation 
in bulk GN gauge. Starting from (3.8), considering the fiy and yy Einstein equations 
(3.14), (3.15), and the mode decomposition discussed in the previous sections we can 
write, 

Sg^^{x,y) = ^/^[0^^(f) + Um,{y)h^^^''\^)+ / dm Ura{y)h^;^\x) 

Sgyy = 6gf,y = , (3.58) 

where (f){x,y) = exp{—Hy /2)(f){x) as before, and the 0^\(j) and h terms automatically 
cancel each other in bulk GN gauge. In order to calculate the effective action, it is 
convenient to fix the brane position so that it lies at y = 0, whilst maintaining the bulk 
GN gauge near infinity. This can be done with a carefully chosen y-dependent gauge 
transformation: 

x^-^x^-i\ y^y-e, (3.59) 



where 



1 for y > i? 1 for y > i? 

(3.60) 

where i? > is some arbitrary finite radius, and a is given in Eq. (3.36). The gauge 
transformation (3.60) should be viewed as the limiting form of a smooth interpolating 

"^ "Forking" , or "dowsing" , is a practice which sometimes reveals an occult presence by means of a 
two-pronged fork, whose role in our case is assumed by the second order effective action. 
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family of test functions that continuously vary in the bulk. The new gauge is not Gaus- 
sian Normal everywhere, interpolating instead between a brane-GN gauge at the brane 
positioned at y = 0, and a bulk-GN gauge near infinity. The bulk metric perturbation 
in this gauge is 

6gab -^ Sgab{x, y) + Va6 + Va6 , (3-61) 

where V is the covariant derivative for gab- The brane metric perturbation in turn is 



57, 



2e 



4 



flU 



Ha 



(3.62) 



where / is the gauge transformation of section 3.3, needed to keep the tensor pertur- 
bation transverse-traceless in the presence of matter perturbations. The second-order 
perturbation of the total action is 



S 



Ml 



[ d'x ^^gdg'^'dGab + \ I d^x v^(57^^ {5Q,, - T^,) 



(3.63) 



and to get the AD effective action we should integrate out the bulk, substituting the 
mode expansion for the radial coordinate y, while keeping all the AD, a;-dependent 
modes off-shell. This means that in the explicit evaluation of the terms in the action 



(3.63) we do not require that (D^ + AH^ 



0, or (D2 -m^ - 2/j2)^M = q. In 
fact, once we have used our Ansatz for the perturbations (3.58), which respects the TT 
conditions, as a means for properly identifying the propagating degrees of freedom in 
the theory about the backgrounds (2.7), we can relax these conditions when working 
out the effective action by evaluating (3.63) on (3.58). The TT conditions for Xi^' will 
nevertheless still emerge from the AD field equations obtained by varying the effective 
action, just like in massive f/(l) gauge theory (see also [68]). We stress that we could 
have used the gauge- fixed action from the start, enforcing TT constraints directly in 
the effective action. We won't do so for the sake of simplicity, because the results are 
completely equivalent at the classical level. 

Using (3.61), (3.62) and the Bianchi identity 'S/"'5Gab = 0, we find that to the 
quadratic order in perturbations the action is 



f d^x^/^ I dy aVah'^'dG^^ - 2M^ f d^x^J^ C{x, 0)5^ 



S = -Ml / dS^/^ 



ay I y=0 



-| / d*x^/^ 



/^^^(x,0) + |o,.(^ + /) 
ri a 



(50 ^u - T^u) 



(3.64) 
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From the metric (3.61) and Eqs. (3.13) - (3.15) it then follows that the variations of 
the Einstein tensor obey 



i\/aSG 



u 



fiu '^ma\yj^'-fj.u 



{y)X^j::^^ + / dm «™(i/)X(-) , 



'AH 

2 



(3.65) 



5G 



fMy\y=0 



Ml 

'ami 



mlu^^{0)D'' {h^;2'^ - h^"^^\,) 



+ / dm m\m{0)D^ {h'p^ - h^'^'>%,) 



(3.66) 
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3eH 



(Mr\ 

2 \2MiJ 



'2- ,(0)/i("^) 



mdUmA 



dm m2M„(0)/i(™) 



3H 
2 



-H\l + e){D^ + 4:H^ 



(3.67) 



-(m) 



In these equations we have been using the tensorial operator X^u , defined by 

X(™) = X^,(/i(™)) + im2 (/^t™) - /.Mv) . (3.68) 

where X^j,(/i*^™)) is given by (3.16). We further use the formula for the variation of the 
brane stress-energy, which is 



se^u = Ml 



«„,(0)X(-'')+ / rfm«„(0)X(-) 

J3|f 

+2 {Ml - MleH) [D^D, - V(^' + 3if ')] / 



(3.69) 



A useful identity which follows from Bianchi identities and stress-energy conservation 
D^'X^,{h) = D^'T^, = is D'^ (50^, - T^,) = -2Ml6G^y\y=o. Using it, the orthogo- 
nality of mode functions Um, the projections 



{ml-2H')(u^^ 



(m^-2H^)(u 



H 

"exp ( -—y 

H 

"exp ( -—y 



-AiMi - MleH)umM 
-A{Ml - MleH)ur,,{0) , 



(3.70) 



and the defining equation (3.49) of the gauge parameter /, after a straightforward albeit 
tedious calculation we finally determine the AD effective action, S'gg = / d'^x^—^ -^eff ' 
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where the Lagrangian density is 



3H_ 

2 



3(1 + e; 3 2 r/n2 , Au2 



M^H' -{D' + AH^)(t) , (3.71) 

and where we use the gauge invariant brane stress-energy perturbation 

r^u = T^. - 4 [Ml - MJeH) [D^D, - %,{D^ + 3H^)] f . (3.72) 

Varying this action reproduces the correct field equations for hf^ and 0. The scalar 
field does not have direct matter couplings at the level of quadratic action, and in 
fact drops out altogether on the normal branch (e = —1), refiecting the fact that the 
normalizable scalar mode in this case is pure gauge. On the self-accelerating branch 
(e = +1), this mode is gauge-invariant, and since it sees the metric it should couple at 
least to gravity at higher orders in perturbative expansion. This mode is itself a ghost 
when a > 0, which occurs when brane tension is negative, as can be readily seen by 
using (2.10). The tensors, /i}I^ , are described mode-by-mode by the standard Pauli- 
Fierz Lagrangian for massive gravity. They couple to matter with the coupling strength 
given by the bulk wave function overlap with the brane 'Um(O). For the continuum 

_3 

modes, this coupling is of the order of M^ ^ . For the discrete mode, it is 



UmM = YT 



1 

M4 



3MlH - 2M|(1 + 
3MiH - 2M|e 



(3.73) 



On the normal branch, this coupling vanishes as o" ^ 0. This is simply the consequence 
that the normalizable zero mode on the normal branch decouples in the limit of infinite 
bulk volume, as is well known [19, 27]. 

In contrast, on the self accelerating branch, the coupling does not vanish, but 
remains of the order of I/M4. As we have already discussed above, for all positive 
values of the tension, the helicity-0 component of this mode is a ghost, which therefore 
remains coupled to matter on the brane with the gravitational strength. It does not 
decouple even when in the vanishing tension limit where the accidental symmetry of 
[43] for the tensor of mass m^ = 2H^ appears, because the symmetry is now realized 
as a Stiickelberg symmetry which mixes the lightest tensor and the normalizable scalar 
mode because they are degenerate. This has also been discussed in the recent work 
[68] . Thus the dynamical degrees of freedom of the tensionless self-accelerating solution 
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are given by the combination (3.45), which we repeat here for completeness, 



where 

The AD tensor H^^ obeys {D^ - AH^)n^,u = -HO^^F. The mixing of the hghtest 
tensor with the brane Goldstone F in effect really just promotes the particular, non- 
dynamical, gauge function / that enforces the TT gauge conditions into a full-fledged 
dynamical mode C, = F + f which mixes with the field Tif^,^. Indeed, in this case we can 
rewrite the boundary condition for 7i^u{x) by subtracting the degenerate eigenmode 
O^pF to (3.48), which will modify the stress-energy source for this mode to 

r^""'^ = T,, + 4M| [D,D, - V(^' + 3//')] i , (3.74) 

where we have used that in the limit of vanishing tension on the self-accelerating branch, 
H -^ 2MI/M'1 and m^ -^ 2H^. Substituting this in the action, and renormalizing 
Ti^u -^ J Ti^y we finally find the zero tension effective lagrangian 

^eff = -\n^''x,,{n) - \H^ {m^n,. - n') + -i T^^-ri^) 



-|- / dm 



3H_ 

2 



_i /.M/^-xi™) + iM™(0)/i(™)'^'^r^J , (3.75) 



where t^^, is given by (3.72), and Tfj,u by (3.74). In this case the scalar ^ does not 
have a normal kinetic term, and only enters the action through mixing with the discrete 
mode TC^u- This occurs because 1/a vanishes as a ^ 0, as per Eq. (3.36). Yet this is 
sufficient to ensure that the ghost survives the vanishing tension limit from the point 
of view of the effective action. 

3.5 Perturbations that break Z2 symmetry 

So far we have been mostly concerned with Z2 symmetric perturbations about the Z2 
symmetric background (2.7), (2.9), following the conventional analysis of the stability 
of DGP backgrounds. However, even with a Z2 symmetric background, if we regard the 
brane as a domain wall rather than an orbifold, which may be well-motivated for DGP 
branes, there is no reason why perturbations about that background should respect 
the Z2 symmetry. Indeed, given the notion that induced curvature is a finite width 
correction for domain walls, one may argue that non-Z2 perturbations are in principle 
just as important as their Z2 symmetric counterparts. 
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To extend the perturbations to non-Z2 symmetric configurations, we take the warp 
factor to still be 

a{y) = exp{eH\y\) , (3.76) 

but imagine that the bulk spacetime describes two separate half-intervals, parameter- 
ized explicitly by — oo < y < and < y < oo, with the brane positioned at y = 0. 
We now consider non-Z2 symmetric perturbations about this background for the case 
where T^^, = 0. These perturbations must satisfy the bulk equations of motion (2.2) 
and the vlovl-TLi symmetric Israel equations (2.3) to linear order. In addition we demand 
continuity of the metric across the brane. This comes for free by construction for Z2 
symmetric perturbations, but not otherwise. If we choose a GN gauge i^^Qyy = Sg^y = 0) 
with brane fixed aX y = 0, we note that the following perturbation satisfies (2.2) and 
(2.3), to linear order, and is continuous at y = 0, 

Sgf^u{x,y) = — — — Va{y) sinh f — J Of,u(j){x) + ^/a{y) / dm sm{ij„,y)x'p\x) 

(3.77) 
where XtJ^ is transverse-traceless, and 

{D' - 2H')x^;:^ = m\^;:\ (D' + AH')^ = (3.78) 

This solution is clearly non-Z2 symmetric since Sgab{x,—y) = —6gab{x,y). Note that 
the tachyonic scalar, 0, is not present for the normal branch (e = —1), but is present 
on the self-accelerating branch (e = +1). It represents yet another instability for the 
self accelerating background. Indeed, for e = +1, (p corresponds to a crinkling up of the 
brane. We can see this by transforming to bulk-GN gauge, which is the appropriate 
gauge for an observer at infinity. To transform to this gauge, we take 

x^^x^ + —^sgn{y), y-^y+- (3.79) 

Ha a 

The bulk metric is now given by 6gyy = 6g^y = 0, and 



^QiJivix.y) = -— sgn(|/)O^^0(x) + \/a{y) / dm sin{ujmy)x^-^\x) (3.80) 

J3fl. 



Since the brane resides y = (f){x) we can immediately interpret the tachyon as a 
crinkling up of the brane. However, we stress that this instability is phenomenologically 
very mild when compared to the ghost, since it is controlled by the time scale given by 
the inverse mass of 0, Tinstabuuy ^ H'^. 
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4. Shocking 4D nonlocalities 

Up until now we have only considered perturbations of the backgrounds (2.7) which are 
normalizable in the bulk. They admit to an effective 4D description mode-by-mode, 
insofar as one is interested in computing their couplings and propagators as measured 
by brane-localized processes. This does not imply that the full picture is 4D over all 
the relevant scales in the infra-red. The additional helicities of massive gravitons spoil 
the 4D description, albeit at very long distances oc r^ by altering momentum transfer 
at very low momenta. As obscure as this may seem in the 4D effective action, it 
becomes transparent in the shock wave analysis of [33]. On the other hand, once one 
restricts only to the normalizable modes localized to the brane, as we have seen above 
one inevitably encounters the lightest graviton on the self-accelerating branch, with 
mass in the unitarity violating window of [44, 45, 46], and so with a helicity-0 ghost. 
This signals an instability which renders the perturbation theory within the effective 
4D description essentially meaningless. Indeed one does not know how to define the 
ground state of the theory on top of which to perturb, and has no clear description 
of the evolutionary end points to which the perturbative ghost may lead. Thus before 
trusting 4D perturbative description one must find ways to neutralize the ghost. 

Since the ghost comes on board with the localized massive graviton, one might try 
changing perturbative definition of the theory, for example by changing prescriptions 
for boundary conditions at infinity, to avoid this mode^. This possibility seems natural 
since after all DGP is really a higher-dimensional theory, disguising as AD at best over 
a finite range of scales. Its ghost arises only after one 'reduces' the theory on the bulk- 
brane background and restricts to the normalizable bulk modes, and so the ghost may 
represent merely an intrinsic instability of this reduction. A brane laden with matter 
may want to move around the bulk in ways which require reintroduction of genuine 
bulk modes, that are not normalizable and hence remain completely outside of the 
scope of the usual 4D effective action analysis. This is supported by the observation of 
[33] that the singular behavior of a shock wave sourced by a photon on the tensionless 
brane may be smoothed by reintroducing a genuine bulk mode, which then resonates 
with the brane. A brane carrying a photon pulse behaves like an antenna, emitting 
bulk gravitons. 

However in this approach, energy will leak from the brane into the bulk as time 
goes on, revealing the fifth dimension. This leakage may eventually strongly deform 
the bulk far away. Alternatively, one could imagine a time-reversal of this process, 
where the description of a particle moving on the brane requires an incoming wave 
in the bulk, with the phase precisely tuned to cancel the ghost divergence, pouring 



^We thank G. Gabadadze for useful discussions about this approach. 
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energy into the brane. One might try cutting the bulk infinity out of the picture, 
seeking boundary conditions that ought to keep the brane stable and self-accelerating. 
It however remains difficult to imagine how this could ever insulate the brane physics 
from distant bulk in the full nonlinear covariant theory, and simultaneously retain the 
guise of a local, causal, 4D description. Every time something happens on the brane, 
one would think that one needs to change the bulk far away, which either requires 
unacceptable external interference, appearing as nonlocalities, or the cross-bulk causal 
transfer of signals that would violate 4D description. Furthermore, while individual 
non-normalizable modes may be treated separately in the linearized theory, they will 
in general couple to each other at higher orders because they will have non-vanishing 
overlaps when nonlinearities are included. Thus once any one non-normalizable mode 
at a fixed 4D mass level is brought back, it should pull alongside it modes at other 
mass levels. These modes may introduce new dangers. 

A full bulk perturbation theory including all non-normalizable modes is beyond the 
scope of our work, requiring first setting up the precise perturbative formulation of the 
problem, defining the set of new orthogonal modes et cetera. However to shed some 
light on the problem we can completely circumvent all those issues by going directly to a 
special limit where we can solve field equations exactly. We shall solve exactly the field 
equations describing the gravitational field of a photon moving on the brane, including 
non-normalizable modes. Such methods have been used in the context of braneworlds 
in [33, 69]. The result of this calculation shows that this time the exchange of the new, 
lightest non-normalizable tensor modes also generates repulsive contributions to the 
potential. This points that light non-normalizable tensor modes behave like 4D ghosts, 
since their contribution to the potential is precisely the off-shell scattering amplitude 
for the single-particle exchange between the brane source and a probe, which in the 4D 
language would be the propagator, that would need to have its sign ffipped to account 
for repulsive force. 

To see this, let us revisit the calculation of [33], describing the shocked background 
geometry (2.7) with the metric 

The induced metric on the brane is 

2 Mudv A5{u)^du'^ 1 — H'^uv 2 d^2 

^^ " (1 + H^uvy " (1 + H^uvy ^ ^ 1 + H^uv ' Ip ■ ^'' 

These metrics are obtained from the static patch form of (2.7) written in terms of the 

1 ll-Hr ^^r^( ZT+\ „^J „, _ 1 ll-Hr 



null coordinates u = -^w } , ^^ exp(Ht) and v = -^w | g-^ exp(—Ht). To determine the 
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field of a plioton in de Sitter geometry, for teclinical reasons it is simplest to actually 
consider the case of two photons with the same momentum p which run in the opposite 
direction in the static patch [70, 71, 72]. So as in [33] we add two antipodal photons on 
the brane, moving along the geodesies u = 0, 9 = 0, and -u = 0, ^ = tt, by introducing 
metric discontinuities along the photon worldlines by substituting dv -^ dv — 6{u)^du. 
This yields (4.1), (4.2). To return to a single source we can multiply this solution by 
the step function 0(7r/2 — 6) as in [72]. The photon stress-energy tensor is 

T^ = -a6^, + 2^g,^J6{9) + 6{9 - 7r))6{<P)6{u)6;6: , (4.3) 
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where we use the notation g4uv for the metric on the brane in Eq. (4.2). A straightfor- 
ward calculation [33] then yields the wave profile equation 



Ml /o2.. , o r.o .. , rr2^. ;. , ...^^ ,...., ...^ .. ^ 2p 



(9j$ + 3e/J9|,|<l> + ij2(A2$ + 2$)) + (A2$ + 2$)5(y) = ^(^6{Q) + 6{Q')y{y) , 

(4.4) 
where we use the shorthand 6{Q) = 6{cos9 — 1)5(0) and 6{Q') = 6{cos9 + 1)5(0). 
The operator A2 is the Laplacian on the transverse 2-sphere on the brane. Using the 
spherical symmetry of the brane geometry transverse to the photon directions, the 
addition theorem for spherical harmonics and linearity of (4.4), we can decompose the 

solution as 

00 

$ = I](*rHl/)^Kcos^) + ^-\y)Pi{- cos 9)^ . (4.5) 

1=0 

Here $1 {z) are the bulk wave functions; $| is sourced by the photon at ^ = 
and $1 by the photon at ^ = vr. By orthogonality and completeness of Legendre 
polynomials, the field equation (4.4) yields an identical differential equation for both 
modes $} (z): 

dl^, + ?>eHd\y\^, + H\2-l{l + l))^, = ^^[^^^:^^ (4.6) 

Using pillbox integration and recalling Z2 symmetry which imposes $/(—!/) = ^i{y) we 
finally determine the boundary value problem for $;: 

dl^i + ?>eHdy^i + H\2- l{l + 1))$, = , 

^i{-y) = <^i{y), (4.7) 

^,, , 2-/(/ + l) , , H21 + 1 
g g 27rM| 
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where g = 2M|/(M|i/) = l/{Hrc). Hence $5+^ = $J ^ = $, since both satisfy the 
same boundary value problem. Further because Pi{—x) = {—iyPi{x), in the expansion 
(4.5) only even-indexed terms survive. This yields 

oo 

<l> = 2Y,'^2i{y)P2i{cose), (4.8) 

1=0 

circumventing the unphysical 4D singularities of / = 1 terms [70, 72]. Solving the dif- 
ferential equation in (4.7) we see that the modes are of the form x ~ e'^^'^'^'^^-'^+^^^/^l^l^L 
In [33] only normalizable mode, for which ||$|p oc f dye^^'^' |$p was finite, were kept. 
Here we want to see what happens when the non-normalizable modes are retained 
instead, and so we use the general bulk wave function 

<^2i = A2,e-(^'+^)^l^l + 52.e(2'-^)^l^l , (4.9) 

where A2i-raode is normalizable and i?2/-mode is not. Substituting this into (4.7), (4.8), 
and introducing the parameter a2i by B21 = —l^TFI7^^^~^ — ^-f'tn , ^ — T+r^0i2h because it 
makes the representation (4.10) particularly transparent, after simple algebra we obtain 



00 _. 
^(n,y) = ^Y(4l + l)P2i(cose)(- TT^^ -, re- 



021 



(21+2) H\y\ 



. ^ A2l-l)H\y\ 



(4.10) 

The parameters a2i are selected by the boundary conditions at the bulk infinity, and, 
in the language of 4D theory, they correspond to the choice of the vacuum, since 
their specification picks out a specific linear combination of the solutions to represent 
a particle state with a given mass and 4-momentum. Clearly, a2i = corresponds 
to keeping only the normalizable modes in the description, and (4.10) reduces to the 
shock wave solution of [33], whereas a2i = 1 selects only the non-normalizable modes, 
throwing out the normalizable ones. 

Now, how should we read the solution (4.10)? First notice that using the spherical 
harmonics addition theorem, ^^^P„(cos^) = X]m=-n^nm(0) 0) ^m(^) 0); ^^^ setting 
n = 21, we can rewrite (4.10) as 

^ OD 21 



4 /=o m=-2l 



X I ^-0^21 g-(2«+2) H\y\ 



(2;-i + ^g)(/ + i + ^g) 

+ (2/-l-i^g)(/ + l-i±^g)" )-^^-''^ 
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This is the Green's function of the problem (4.4), describing the gravitational field of 
a 'particle' of effective mass p in the space transverse to the photon's u, v propagation 
plane. Now, we recall that in the conventional approach, a tree-level perturbative 
potential generated by an exchange of a mediating boson is the Fourier transform 
of the scattering amplitude involving the boson propagator. The formula (4.11) is 
precisely the same: since the transverse space to the shock is a S^^^^^ x Rbuik, the 
terms of the expansion (4.11) can be interpreted as the off-shell tree level amplitudes 
involving the exchange of the discretized Fourier modes on the sphere, described by 
the spherical harmonics, and weighed by the bulk radial functions which account for 
the dilution of intermediary's wave function due to the warping of the bulk, when the 
intermediary slides just outside of the brane. The 'quantum numbers' / measure the 
Euclidean momentum of the intermediary modes on the 3"^^^^^^, q r^ HI, and control 
the momentum transfer in a scattering process between the photon and a distant test 
particle mediated by a virtual intermediary. The bulk y-dependence scales the coupling 
up or down depending on the location of the probe versus the brane. 

Now, to see how the shock wave looks on the self-accelerating brane, we can set 
e = 1 and y = in Eq. (4.10), which yields 

„ oo 2/ _. 

(4.12) 
These formulae are very revealing. First, let us suppose that < 0:2; < 1. It is clear 
from the y-dependence of (4.10) that general solutions with a2i 7^ peak far from the 
brane, indicating that they are very sensitive to the perturbations near the bulk infinity. 
However, from (4.12) we see that the local physics on the brane is not that sensitive to 
distant bulk, since for large momentum transfer / ^ g, i.e. at short distances between 
the source and the probe, HTZ ~ a/2(1 — cos^) <^ 1, the a2i terms which encode bulk 
boundary conditions completely cancel in the leading order in (4.10), (4.11), (4.12): 

( (2l3iTiKmy + i2i-i)T+i-i) ) = ^ + ^(f)- '^^^^ *° *^^ leading order (4.12) behaves the 
same at short distances as the wave profile with 0:2/ = 0, composed only of normalizable 
modes. Its short distance form will be very well approximated by the Aichelburg-Sexl 
wave profile, in exactly the same way as the purely normalizable wave profile [33]. 

However, at very low momentum transfer, or at very large distances, a general 
solution (4.10) with non- normalizable modes differs very dramatically from the one 
built purely out of the normalizable modes. The point is that the coefficients of the 
expansion of the wave profile oc a2i are not positive definite when viewed as a function 
of /. Indeed, when g > 2, all the coefficients in the second term change sign for 
all / < g/2 — 1 except for / = 0, which remains positive because of the 1/(2/ — 1) 
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factor. Now when g < 2, terms with / > 1 remain positive, but the / = term turns 
negative. Hence these modes come in with opposite signs relative to the contributions 
from their normahzable partners oc 1 — a2i, and because their dependence on the 
transverse distance from the source is the same as for the normahzable modes, this 
means that they yield repulsive contributions to the gravitational potential at large 
distances. Given our interpretation of the contributions of the terms in the expansion of 
the wave profile (4.11) as the discretized propagator of the exchanged virtual graviton, 
the repulsive contributions to the potential at very large distances signal a ghost-like 
behavior among the lightest states in the non-normalizable tensor sector. Note further 
that when g is an integer, there are pole-like singularities in the sums (4.10)-(4.12). In 
the normahzable sector, the only pole in fact occurs when g = 1, that corresponds to 
the limit of vanishing tension and as discussed in [33] which can be interpreted as the 
instability that leads to rapid release of energy into the bulk, related to the lightest 
tensor ghost. The non-normalizable contributions however have poles at all even integer 
values of g > 0, that indicate that if energy were to be inserted into them at infinity, it 
would be quickly transferred into the normahzable modes, which could produce large 
gravitational effects on the brane. 

This shows that resurrecting non-normalizable modes, if viewed as 4D phenomena, 
may open the door to new ghost-like modes, over and above the helicity-0/scalar ghost. 
In this case one further needs to recheck carefully the scalar sector for any pathologies 
among the non-normalizable scalar modes, which cannot be revealed by the shock 
wave analysis as they are not sourced by relativistic particles. Note that taking a2i 
outside of the interval [0, 1] will not remove the repulsive contributions to the shock 
wave, but would further exacerbate the problem. For a2i > 1, all the normahzable 
mode contributions would switch sign, whereas for a2i < 0, all but the lightest non- 
normalizable modes would be repulsive. Thus it appears that the repulsive terms in 
the non-normalizable tensor sector will be avoided only if we set a2i = 0, but then this 
retains only the normahzable modes in the helicity-0 sector as well, leaving one with 
the helicity-0 ghost. 

Interestingly, one can check explicitly using (4.10) that the situation on the normal 
branch is better, in the sense that the relative sign between the normahzable and non- 
normalizable contributions remains the same for all the terms in the expansion. In 
fact, on the normal branch, the only term which contributes to the wave profile with 
a negative sign is the / = mode, but this mode is completely constant on the brane 
and produces no force on brane particles. It will only repel bulk probes, but this is 
consistent with the picture that domain walls in flat space, when viewed from the wall's 
rest frame, exert repulsive force on particles in the bulk [55]. 
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5. Summary 

In this work we have reconsidered the perturbative description of codimension-1 DGP 
vacua. Our results confirm that in the normahzable sector of modes, there is always 
a perturbative ghost on the self-accelerating branch. For positive brane tension, it 
resides in the localized lightest graviton multiplet as the helicity-0 state, whereas for 
the negative tension it is the scalar 'radion'-like field. In the borderline case of zero 
tension, describing a background where the brane expansion accelerates solely due to 
gravity modification, the ghost is an admixture of the helicity-0 and scalar modes, 
which become completely degenerate. In contrast, the normal branch vacua are free of 
ghosts, at least when the brane tension is non-negative. 

The ghost makes a simple 4D perturbative analysis of the self-accelerating dy- 
namics prohibitive. It signals an instability which renders the effective 4D description 
perturbatively meaningless. Off hand, one does not know how to define the ground 
state of the theory, and has no clear description of the evolutionary end points to 
which the perturbative ghost may lead. To illustrate the dangers from ghosts, let us 
review a very simple, intuitive example of the classical ghost instability which, to our 
surprise, does not seem to be widely discussed. Consider a system of two degenerate 
harmonic oscillators x and y. One could solve it by solving each individual oscillator 
problem and then define the full set of states as the direct product of the individual 
oscillators. However, these states do not naturally refiect the 0(2) rotational symmetry 
of the system, which can be made manifest by going to polar coordinates x = TZcosO, 
y = T^sin^. Then the angle ^ is a Goldstone mode, whose conserved charge is the 
angular momentum of the system, while the radial variable is an anharmonic oscilla- 
tor moving in an effective potential composed of the original parabolic piece far away 
and the centrifugal barrier erected by the angular momentum close in. This system is 
classically stable. Now imagine making one of the original harmonic oscillators purely 
imaginary, e.g. y -^ iy. This changes the symmetry group to 0(1, 1), the 'polar' coor- 
dinate maps become hyperbolic functions, and the phase is now the ghost, whereas the 
radial mode remains a normal field. The ghost's centrifugal contribution to the radial 
motion is now an infinite well, instead of a barrier. Thus even a tiniest perturbation 
with non-zero angular momentum will send the oscillator spiralling into the infinitely 
deep centrifugal well, spinning it up indefinitely as it goes in. In terms of the origi- 
nal Cartesian variables, the two oscillators gain infinite kinetic energy at the expense 
of each other, while keeping their difference constant. When the physical oscillator x 
couples to other normal degrees of freedom, it can transfer its kinetic energy to them, 
destabilizing the rest of the world. The rate of the instability is controlled by the os- 
cillator period, and in a field theory where one has a tower of ghastly oscillators with 
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arbitrarily high frequencies, the energy transfer rates can therefore be very fast. Inter- 
estingly, this situation is reminiscent^ of what one encounters on the self-accelerating 
branch of DGP, where among the degenerate lightest gravitons there is a ghost. 

We should stress that we do not claim that it is ultimately impossible to get rid 
of the ghost. The question is, what is the price one must pay to exorcize it. We have 
considered what happens when one restores non-normalizable modes in the bulk, and 
allows them to couple to brane matter. This may be an interesting arena to explore if 
by abandoning description in terms of normalizable modes alone, some of the effects 
of the helicity-0 ghost may be controlled. If one reinstates the bulk scalar modes that 
could couple to the helicity-0 ghost, however, by bulk general covariance one should 
also bring in the tensors. All of these modes dwell outside of the realm of 4D effective 
theory at all scales. Even so, our direct calculation of their couplings to relativistic 
brane matter, modelled by a photon zipping along the brane, shows that at short 
distances from it, the photon's gravitational field follows closely the 4D form, being 
indistinguishable from the purely normalizable contributions in the leading order of the 
expansion in g = I/Htc- However, we have noted new repulsive gravitational effects at 
large distances from the source, arising from the non-normalizable tensor sector, which 
indicate that low momentum tensors behave as ghosts. This issue clearly deserves more 
attention, and it would be interesting to develop a description of the processes beyond 
quadratic perturbation theory to see how the system evolves. 

We note in passing that another possibility may be to consider altering the theory 
at the level of the action itself, for example by adding extra terms in the bulk or on the 
brane, and compactifying the bulk by adding new branes^°. In the former case adding 
the Gauss-Bonnet term in the bulk seems an intriguing possibility since the resulting 
field equations are of second order and so they maintain the simple distributional 
brane setup without adding new, possibly dangerous, graviton states of generic higher- 
derivative models. Cutting the bulk at a finite distance will reintroduce the zero mode 
graviton, and may change the boundary conditions for the massive modes, affecting 
the values of masses. If this could lift the localized lightest tensor out of the unitarity- 
violating window, and/or break degeneracies with scalar modes, the ghost instability 
may be tamed. One then expects to be left with a strongly coupled massive graviton 
on the original DGP brane, which could dominate over the zero mode in a range 
of scales, and so one may consider phenomenological aspects of such a multi-gravity 
theory. However compactifying the bulk would turn DGP immediately into an effective 



^The differences are tlie presence of background de Sitter geometry and a different isospin group 
(0(2, 1) instead of 0(1, 1)) but much of the rest appears the same, at least the hncarizcd perturbation 
theory level. 

"'^'^We understand that K. Izumi and T. Tanaka are pursuing such approaches. 
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AD theory at large scales, exposing it to the edge of the Weinberg's venerable no-go 
theorem [4] for the adjustment of cosmo logical constant. Thus while removing the 
ghost by a compactification of the bulk might work, it may automatically restore the 
usual fine-tunings of the AD vacuum energy, completely obscuring the whole point of 
self-acceleration. 

Thus it is doubtful that self-acceleration in its present guise may serve as a model 
of the current epoch of cosmic acceleration, since after all it does appear that some 
perturbative description of our universe at the largest scales should exist. The self- 
accelerating branch does not seem to fit this bill due to its occult sector. Yet, if 
one wants to study the implications of modified gravity, one may still find a useful 
framework among the brane-induced gravity models. The simplest one may be the 
normal branch solutions. True, they undergo cosmic acceleration at the right rate 
because one fine-tunes the brane tension to just the right value by hand, a ~ (lO'^eV)'^. 
But there is no ghost, and perturbative description is reliable. If one then also tunes the 
scale of modification of gravity, Vc ~ ^/H, one gets interesting signatures of weakened 
gravity at the largest scales. Namely although there is a zero mode graviton on the 
normal branch, the light bulk modes also contribute to gravity at scales smaller than 
Tc ~ ^/H and the effective graviton that mediates sub-horizon interactions is really a 
resonance composed of many modes. Thus local gravity is stronger than the horizon 
scale gravity. At the horizon scale, the force weakens because the effective momentum 
transfer of the massive admixtures changes from 1/g^ to 1/q as the extra dimension 
opens up. This weakening of gravity may simultaneously change the cosmic large scale 
structure [34] and masquerade the cosmo logical constant as dark energy with w < —1 
[73], in a way that could be accessible to observations. 
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6. Proof of tensor decomposition in Eq. (3.8) 

To prove (3.8) we first introduce auxiliary fields S^^, = h^y — j'j^u, B^ and _B, where 
we require that B^ and B are solutions of differential equations 

{D^ + AH'')B = '^D'^D'^^, , 

D^B^ + D^'D^B, - ^D^B = D^S^, . (6.1) 

These differential equations are integrable, despite the cumbersome nature of the 
vector field equation that appears to mix different components. We can simplify the 
system by introducing an additional scalar auxiliary field, as follows. First, commute 
through the derivatives acting on the vector, using the standard rules for commutators 
of covariant derivatives in de Sitter metric 7^jy and rewrite the vector equation as 

{D^ + 3H^)B^ = D"^^, - D^iD.B" - ^B) . (6.2) 

Next multiply this equation by D^, commute the derivatives again, and using the 
defining equation for B to eliminate D^D^H^y, finally obtain the equation for \E' = 
B - D^B^' : 

{D^ + 3H^)^ = 0. (6.3) 
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Thus \1/ is a completely free field on de Sitter background, and this equation can be 
solved at least in principle. 

We can now consider a different system of equations governing the auxiliary fields, 
where we eliminate D^B^ in (6.2) in terms of ^ and B. The full system then becomes 



(Z}2 + 3ij2)^ __ 


= 0, 






{D^ + 4:H^)B -- 


= -Df^D" 


'y 




{D' + 3H^)B^ -- 


= D'l:^, 


-D,{ 


B 

~2 



^). (6.4) 

Clearly, once we find a solution \1' and determine a B, which we can do because we know 
the source in the B equation, given by the original tensor h^^, we can integrate the 
remaining 4 equations for the vector to get the full solution. It will, clearly, depend on 
the choice of the auxiliary function \E'. To select precisely the required solution of (6.1) 
we must extract that solution of (6.4) which automatically satisfies ^ = B — Df^B^^. So 
indeed multiplying the vector field with V^ yields 

3i7=^^ = (D^ + QH'^) {B - D^B^) , (6.5) 

and when we exclude the null eigenvalue spurious of D^ + QH^ from B — D^B^ by 
an appropriate choice of boundary conditions'^, the field equation {D^ + 3i/^)^ = 
implies (D^ + 3H'^){B - D^B^") = 0. In this case ^^^ = B - D^B'' solves the same 
equation as ^. But then, since (D^ + 6if^)\E'i = 3iJ^^i, comparing with (6.5) gives 
\l'i = ^. Thus as required "$ = B — D^B^, and any such solution of (6.4) will be 
exactly a solution of (6.1). Having a solution, we can now define the tensor 

hf,u = S^j, - D^B,, - D^B^ + — %^, , (6.6) 

and note, using the differential equations (6.1) that it is transverse, Dfji^y = 0. Its 
trace is 

/i% = 2B- 2D^B>' , (6.7) 

and so we can eliminate B from the solution writing it as -B = D^^B^ + h/2. Moreover, 
we can separate the vector field B^ as a Lorentz-gauge vector A^, D^A^ = 0, plus a 
scalar gradient. 



B^ = A, + -D^(f), {6.1 



^^We can always exclude these spurions, because by {D^ + 01^)0 ^d = D^{D'^ + QH'^)d, the gauge 
shifts B^ — > i?^ + D^§ by functions obeying {D^ + QH'^)d — drop out from the field equations (6.4) 
but shift the spurion in 5 — D^B^ by 6-ff^i?. So we can simply choose ^ to completely cancel away 
the spurion. 
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and substitute all this back in (6.6). Solving for the original tensor field /i^jy, we write 

I h — h 

V = V + D^,A^ + D^Af, + D^D^(f) - - V^ V + ^^ V > (6-9) 

Now, we are almost there: the field /i^jy is transverse, but not yet traceless. However 
its trace is equal to the auxiliary free field ^ by Eq. (6.7), which we can pick to be 
exactly zero by choosing appropriate boundary conditions, amounting to gauge fixing 
B^. Hence the tensor is 

Vk=o = ^"> (6-10) 

which is also traceless! Thus indeed we recover (3.8), as claimed. Note that nowhere 
in this decomposition did we need to specify anything about y-dependence, which 
was treated as an extra parameter. Hence, Eq. (3.8) will remain valid for a Fourier 
transform of the metric perturbation as well. So we see that an arbitrary perturbation 
in the GN gauge can be separated as a 5-component transverse-traceless tensor, a 3- 
component Lorentz vector, and 2 scalars, in addition to the brane location F, which is 
a separate AD field. 

7. Helicity-0 ghosts in de Sitter space 

Here we review the proof that a helicity-0 mode of the massive Pauli-Fierz spin-2 
theory in de Sitter space is a ghost if the mass obeys < m^ < 2H^. In the early 
work of Higuchi [44] this was demonstrated by showing that the helicity-0 sector of the 
Hilbert space contains negative norm states, but since then simpler methods based on 
Hamiltonian analysis have been developed [46]. Here we largely follow the analysis of 
[46], although we note that a Lagrangian analysis based on correctly identifying the 
residual gauge symmetries of the helicity-0 sector would produce equivalent results. 

We start with the Pauli-Fierz massive spin-2 theory in a background metric 7^j,, 
which is given by £ = v^— 7 — \h^^X)j^\ where X)j^ is defined in Eqs. (3.68) and 
(3.16). We take ^^^ to be the de Sitter metric (2.8), and D^ its covariant derivative. 
This Lagrangian describes the localized lightest tensor multiplet on the self-accelerating 
branch, with mass m = rrid, which for positive brane tension lies in the unitarity- 
violating window < m^ < 2H^. The explicit form of the Pauli-Fierz action is 

1 



-iH^h^'^h,, - iH^h' - ^ {h^^-h,, - h') , (7.1) 

where the /i^jy are the general metric perturbations. Now, using the after-the-fact 
wisdom [46], we know that if the mass m'^ were zero, this theory would only have two 
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helicity-2 excitations as the propagating modes. Thus the modifications can only arise 
because some of the scalar and vector perturbations of the general 4D metric do not 
decouple when rn^ ^ 0. Further, because the theory remains Lorentz-invariant, the 
vectors and the scalars decouple from the tensors, and moreover the vectors can only 
yield the helicity-1 modes. Thus the helicity-0 mode can only arise from the scalar 
perturbations, which we can parameterize as 

h,, = 2e-5^* {d,d,E + A6,,) , hu = e-^^* d,B , hu = 26"!^* . (7.2) 

Here di are spatial derivatives. We have normalized the perturbations by the appropri- 
ate powers of c^* to simplify the analysis, following [46]. Plugging the perturbations 
(7.2) into (7.1) yields, with Spf = J dtd^x £, 

C = -6A^ - aAx + (4m2 - 9H^)AX + A Um^ - —H^ - 26-^^* A j A 

AHX + 2(m2 - 3iJ^)X + 12HA + Q{m^ - 3H^)A - Ae~^"^AA 
+e-"^AB (aH(J) + aA- 6Ha) - QH^(j)^ - hn^BAB , (7.3) 



where we have introduced X = AE, with A = Yli df the fiat 3-space Laplacian, and 
denoted time differentiation by a dot. To define the Hamiltonian as in [46], we write 
down the conjugate momenta IIq = |^, which are 

Ux = -AA - AH(j) , n^ = -12i -AX- 12if0 + Ae'^^^AB , 

Ub = U^ = 0. (7.4) 

Thus the Lagrangian and the Hamiltonian are related by £ = HxX + HaA — 7i, and 
so using this and (7.4) the Hamiltonian is given by 

n = -UJ,- -UxUa + A (—H^ - 6m' + 2e-''''A^ A + {9H' - Am') AX 

+0 [2{m' - 3H')X + 6{m' - 3H')A - Ae'^^'AA - HUa] 

+e-"^AB {lix + QHA) - \m'BAB (7.5) 

From the Hamiltonian we immediately see that the fields and B aren't propagating - 
which of course comes as no surprise, since they are the scalar remnants of the shift and 
lapse functions. Varying the Hamiltonian with respect to them yields the Hamiltonian 
and momentum constraints respectively, which can be written as 

Ha - — [2{m' - 3H')X + 6{m' - 3H')A - Ae^'^'AA] = , 

-Ht 

B + ^{Ux + QHA) = , (7.6) 

m,^ 
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where we are holding m^ 7^ in the last equation. Substituting these equations in the 
Lagrangian we integrate out the lapse and shift, and integrating by parts find 



c^,n.-li,„'-3HM)x-uJl-i^_ 



2HtA\ 1 

Ux + —im^-3H'')UxX 



+n 



X 



^{m^-3H^)-^{m^-6H^)e-^-^A 



A + (4m2 - 9H'^)AX 



+A 



Qm 



27 



H' 



6 



m^ 



(m^ - 3H^)e 



-2Ht 



A 



A 



(7.7) 



A field redefinition Ilx = p + -jj(jn^ — 3H^)A recasts the Lagrangian as 



C = pX — p 



3 e-2^*A\ e-2^*A . z/^ SmV^ ^^,,2 



2171^ 



p + p- 



H -^^m^'' 



2^.^^ + ^A-^. P^* 



where v^ 



m 



2H^. This shows that with these variables A is not a dynamical 



field. Its field equation is algebraic, and for m? ^ 2H'^ it yields A = —-^ ~ "^^i^i" P- 
Substituting this into (7.8) gives 



C = pX — p 



■^~4HtQA g-2m^ 3 



6m? u"^ 



2171^ 



p — p 






3H 



2H 



2,, 2 



X 



m V 



6iJ2 



X\ 



At long last, we make the one last field redefinition 

H 



X = q + 



2rn?v'^ 
which casts the Lagrangian in the form 



[3(m2-3if2)_2e-2^*A]p, 



C = pq — 



■q - 



6/72 



2m? v"^ 



9H^ 



p m 



.-2Ht 



A p. 



(7.9) 
(7.10) 

(7.11) 



This equation looks slightly unusual, since the Lagrangian seems to depend on the 
spatial gradients of the 'momentum' p rather than the 'field' q. However, this is just a 
mirage, which can be easily removed by a canonical transformation q = — tt, p = ip, and 
the integration by parts of pq = —ipir = nip — -^(ipn). Dropping the total derivative, 
we can extract the final Hamiltonian from the Lagrangian C = Tup — 7i, to find 



2,, 2 



n 



m V 



6H^ 



■Tl' + 



3IP 

2m'^v'^ 



9H^ 



(p \m, — 



— e 



A v^. 



(7.12) 



We immediately see that when u^ < (i.e. < m^ < 2H'^) the Hamiltonian is 
negative definite, and so the field v^ is a ghost. It can be viewed as literally a massive 
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scalar field covariantly coupled to de Sitter gravity, with the sign of the Lagrangian 
reversed. When m^ = and z/^ = the ghost decouples in the pure Pauli-Fierz 
theory, which can be glimpsed at from the canonically normalized scalar cp = ^^^^Lfc, 

n = ^^TTfj indicating that all the perturbative Lagrangian couplings of ip to matter 
are proportional to positive powers of m|z/|. This does not occur for the self- accelerating 
branch DGP because of the additional scalar localized mode, as discussed in the text 
and in [68]. 

Also note that the < rn^ < 2H^ ghost has a very mild tachyonic instability, 
induced by de Sitter expansion. One can see it from the field equation for the Fourier 
components of ipc, which from (7.12) is, by setting ipc = 'fc{k)e^^'^, 

0c{k) +(m^-^ + ee-^'^'jMk) = . (7.13) 

As time goes on, the 3-momentum term oc /c^ becomes insignificant, so that 0cik) -^ 
exp(±-i/^^^ m"^ t) for m^ < 2H^. This behavior changes for all gravitons above de 



Sitter gap m'^ > ^-^ (as discussed in [46]), and when m^ = the ghost is absent in the 
first place. Yet it is clear that these instabilities simply correspond to the freezing out 
of long wavelength ghost modes at super-horizon scales, as is common in inflation, and 
it would be interesting to explore the implications of this mechanism for DGP. 
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